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1 INTRODUCTION: 

In 1953 the physicists E. Inonii and E.P. Wigner (1963 Nobel prize) introduced the concept of 
deformation of a Lie algebra by considering the composition law (u,v) — > (u + v)/(l + (uv/c 2 )) 
for speeds in special relativity (Poincare group) when c is the speed of light, claiming that the 
limit c — > oo or 1/c — > should produce the composition law (u,v) — > u + v used in classical 
mechanics (Galilee group) ([22]). However, this result is not correct indeed as 1/c -> has no 
meaning independently of the choice of length and time units. Hence one has to consider the di- 
mensionless numbers u — u/c, v = v/c in order to get (u, v) — > (u + v)/(l + uv) with no longer any 
pertubation parameter involved ([32]). Nevertheless, this idea brought the birth of the theory of 
deformation of algebraic structures ([13], [19], [20], [21], [34], [35], [44], [45]), culminating in the use of 
the Chevalley-Eilenberg cohomology of Lie algebras ([9]) and one of the first applications of com- 
puter algebra in the seventies because a few counterexamples can only be found for Lie algebras 
of dimension > 11 ([2]). Finally, it must also be noticed that the main idea of general relativity 
is to deform the Minkowski metric dx 2 + dy 2 + dz 2 — c 2 dt 2 of space-time by means of the small 
dimensionless parameter 4>/c 2 where <f> = GM/r is the gravitational potential at a distance r of a 
central attractive mass M with gravitational constant G ([11], [33]). 

Let G be a Lie group with identity e and Lie algebra Q — T e (G), the tangent space of 
G at e. If a = (a T ) with r = l,...,p are local coordinates on G, the bracket [Q,G] C Q is 
defined by p 2 (p — l)/2 structure constants c = (c T pa = —c T ap ) satisfying the Jacobi identities 
J(c) = c^ T + c^ r c^ p + cPpC^ = 0. If now c t = c + tC + ... is a deformation of c sat- 

isfying J{c t ) = 0, the idea is to study the vector space dJ ^ C = modulo the fact that c behaves 
like a 3-tensor under a change of basis of Q. The first condition implies C G Z 2 (G) as a cocycle 
while the second implies C £ B 2 {G) as a coboundary, a result leading to introduce the Chevalley- 
Eilenberg cohomology groups H r (Q) — Z r (Q) / B r (Q) and to study Bi2(G) in particular. 



A few years later, a deformation theory was introduced for structures on manifolds, generally 
represented by fields of geometric objects like tensors and we may quote riemannian, symplectic, 
complex analytic or contact structures with works by H. Goldschmidt ([12], [14], [15]), R.E. Greene 
([15]), V. Guillemin ([16], [17], [18]), K. Kodaira ([24], [25]), M. Kuranishi ([27]), L. Nirenberg ([24]), 
D.C. Spencer ([15], [25], [47], [48]) or S. Sternberg ([16], [17], [18]). The idea is to make the underlying 
geometric objects depending on a parameter while satisfying prescribed integrability conditions. 
The link betwen the two approaches, though often conjectured, has never been exhibited by the 
above authors and our aim is to sketch the solution of this problem that we have presented in 
many books, in particular to study the possibility to use computer algebra for such a purpose. It 
must be noticed that the general concept of natural bundle and geometric object is absent from 
the work of Spencer and coworkers ([27], [47], [48]]) though it has been discovered by E. Vessiot as 
early as in 1903 ([51]). It must also be noticed that the introduction of the book "Lie equations" 
([26]) dealing with the "classical" structures on manifolds has nothing to do with the remaining 
of the book dealing with the nonlinear Spencer sequence and where all the vector bundles and 
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nonlinear operators are different (See ([43]) for more details). This work is a natural continuation 
of symbolic computations done at RWTH- Aachen university by M. Barakat and A. Lorenz in 2008 
([3],[30]). 

The starting point is to refer to the three fundamental theorems of Lie, in particular the third 
which provides a way to realize an analytic Lie group G from the knowledge of its Lie algebra 
Q and the construction of the left or right invariant analytic 1-forms on G called Maurer-Cartan 
forms lj t — (c*£ (a) da' 7 ) satisfying the well known Maurer-Cartan equations duj T + c T pa ui p A u) a = 
where d is the exterior derivative, provided that J(c) = 0. 

The general solution involves tools from differential geometry (Spencer operator, J-cohomology) , 
differential algebra, algebraic geometry, algebraic analysis and homological algebra. However, the 
key argument is to acknowledge the fact that the Vessiot structure equations ([51], 1903, still un- 
known today !) must be used in place of the Cartan structure equations ([4], [5]), along the following 
diagram describing the reference ([43]) which is the latest one existing on this subject while pro- 
viding applications to engineering (continuum mechanics, electromagnetism) and mathematical 
(gauge theory, general relativity) physics. 

CARTAN — > SPENCER 

s 

LIE X X 

\ 

VESSIOT — > JANET 



In a word, one has to replace Lie algebras by Lie algebroids with a bracket now depending 
on the Spencer operator and use the corresponding canonical linear Janet sequence in order to 
induce a new sequence locally described by finite dimensional vector spaces and linear maps, even 
for structures providing infinite dimensional Lie algebroids (contact and unimodular contact struc- 
tures are good examples as we shall see). The cohomology of this sequence only depends on the 
"structure constants " appearing in the Vessiot structure equations (constant riemannian curva- 
ture is an example with only one constant having of course nothing to do with any Lie algebra !). 
Finally, the simplest case of a principal homogeneous space for G (for example G itself as before) 
gives back the Chevalley-Eilenberg cohomology. 

In order to motivate the reader and convince him about the novelty of the underlying concepts, 
we provide in a rather self-contained way and parallel manners the following four striking examples 
which are the best nontrivial ones we know showing how and why computer algebra could be used 
and that will be revisited later on. We invite the reader at this stage to try to imagine any possible 
link that could exist between these examples. 

EXAMPLE 1.1: (Principal homogeneous structure) When T is the Lie group of transformations 
made by the constant translations y l = x l + a 1 for % = 1, n of a manifold X with dim(X) = n, 
the characteristic object invariant by T is a family ui = (uj t = uijdx 1 ) of n 1-forms with det(u>) 
in such a way that F = {/ e aut(X)\ji(f)~ 1 (iv) = uj} when aut(X) denotes the pseudogroup of 
local diffeomorphisms of X, j q (f) denotes the derivatives of / up to order q and ji(f) acts in the 
usual way on covariant tensors. For any vector field £ € T = T(X) the tangent bundle to X, 
introducing the standard Lie derivative £(£) of forms with respect to £, we may consider the n 2 
first order Medolaghi equations: 

fir ee = 0J T r (x)diC +?d r u>T(x) = 

The particular situation is found with the special choice ui = (dx l ) that leads to the involu- 
tive system d£ k = 0. Introducing the inverse matrix a = (a l T ) = uj^ 1 , the above equations 
amount to the bracket relations [£, a T ] — and, using crossed derivatives on the solved form 
di£ k + a^(x)d r uij (x) = 0, we obtain the n 2 (n — l)/2 zero order equations: 
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r0r(^(*)4(aO(0i"JC»O - «(*))) = 
The integrability conditions (IC) are thus the n 2 (n — l)/2 Vessiot structure equations: 

diuj^x) - djujj(x) = c r pa uj?{x)LJ?{x) 

with n 2 (n — l)/2 structure constants c = (c T = —c T ap ). When X = G, these equations can be 
identified with the Maurer-Cartan equations (MC) existing in the theory of Lie groups, on the con- 
dition to change the sign of the structure constants involved because we have [a p ,Q! CT ] = — c pa a T . 
Writing these equations in the form of the exterior system du T = c^ a uj p A u a and closing this 
system by applying once more the exterior derivative d, we obtain the quadratic IC: 

c A c M + c A c M + c x r M =0 

also called Jacobi relations J(c) = 0. Finally, if another family ui of forms is given, the equivalence 
problem = ui cannot be solved even locally if c ^ c. 

EXAMPLE 1.2: (Riemann structure) If ui = (ujij = uiji) is a metric on a manifold X with 
dim(X) = n such that det(uj) ^ 0, the Lie pseudogroup of transformations preserving oj is 
r = {/ G aut(X)\ji(f)~ 1 (lu) = cu} and is a Lie group with a maximum number of n(n + l)/2 
parameters. A special metric could be the Euclidean metric when n = 1, 2, 3 as in elasticity theory 
or the Minkowski metric when n — 4 as in special relativity. The first order Medolaghi equations: 

flij EE (£(£)u})ij EE LO rj {x)di£ r + LO ir {x)dj^ r + C^r^i] (%) = 

are also called Killing equations for historical reasons. The main problem is that this system is 
not involutive unless we prolong it to order two by differentiating once the equations. For such a 
purpose, introducing uj^ 1 — (a/- 7 ) as usual, we may define the Christoffel symbols: 

= ^uj kr (x)(d l uj rj (x) + djU ri {x) - d T u)n(x)) = rfi(x) 

This is a new geometric object of order 2 providing an isomorphism ~ (a;, 7) and allowing 

to obtain the second order Medolaghi equations: 

I* ee (£(07)?j = d l3 e + l k rj {x)diC + li{x)d 3 C - ^(x)d r Z k + CdrlU x ) = 
Surprisingly, the following expression: 

p%{x) ee d a %{x) ~ d n k u {x) + Jlj(x) 7 k ri (x) - J^xh^ix) 

is still a first order geometric object and even a 4-tensor with n 2 (n 2 — 1)/12 independent compo- 
nents satisfying the purely algebraic relations : 

Plij + Pijl + Pjli = 0, UrlPkij + UkrPlij = 

Accordingly, the IC must express that the new first order equations i?^ ■ ee (£(£,)p)iij = are onr y 
linear combinations of the previous ones and we get the Vessiot structure equations with the only 
structure constant c: 

Plij(x) = c(<S?wy(a;) - 5!jwu(x)) 

describing the constant Riemannian curvature condition of Eisenhart ([10]). If Co is another non- 
degenerate metric with structure constant c, the equivalence problem ji(f)~ 1 (co) = Co cannot be 
solved even locally if c ^ c. 

EXAMPLE 1.3: (Contact structure) We only treat the case dim(X) = 3 as the case dim(X) = 
2p+l needs much more work ([37]). Let us consider the so-called contact 1-form a = dx 1 — x 3 dx 2 
and consider the Lie pseudogroup T C aut(X) of (local) transformations preserving a up to a 
function factor, that is T = {/ 6 aut(X)\ji(f)~ 1 (a) — pa} where again j q (f) is a symbolic way 
for writing out the derivatives of / up to order q and a transforms like a 1-covariant tensor. It 
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may be tempting to look for a kind of "object " the invariance of which should characterize T. 
Introducing the exterior derivative da = dx 2 A dx 3 as a 2-form, we obtain the volume 3-form 
a A da = dx 1 A dx 2 A dx 3 . As it is well known that the exterior derivative commutes with any 
diffcomorphism, we obtain sucessively: 

ji(/) _1 (da) =d(ji(/) _1 (a)) = d(pa) = pda + dp A a => (a A da) =p 2 (aAda) 

As the volume 3-form a A da transforms through a division by the Jacobian determinant A = 
9(f 1 ,f 2 ,f 3 )/d(x 1 ,x 2 7 x 3 ) ^ of the transformation y = f(x) with inverse x = ,f^ 1 (y) = g{y) 7 
the desired object is thus no longer a 1-form but a 1-form density u = (uji,lj 2 ,uj 3 ) transforming 
like a 1-form but up to a division by the square root of the Jacobian determinant. It follows that 
the infinitesimal contact transformations are vector fields £ G T = T(X) the tangent bundle of X, 
satisfying the 3 so-called first order Medolaghi equations: 

SU = = Ur{x)9iC - 0-/2)Ui{x)drC + ? d r OJi{x) = 

When ijj = (1, — x 3 , 0), we obtain the special involutive system: 

d 3 e + d 2 e + 2x 3 o 1 e - d,e = o, o 3 e - x 3 o 3 e = o, o 2 e - x 3 o 2 e + x 3 o 1 e - (x^d^ 2 -e = o 

with 2 equations of class 3 and 1 equation of class 2 obtained by exchanging x 1 and x 3 (see section 
4 for the definitions) and thus only 1 compatibility conditions (CC) for the second members. 
For an arbitrary co, we may ask about the differential conditions on co such that all the equations 
of order r + 1 are only obtained by differentiating r times the first order equations, exactly like in 
the special situation just considered where the system is involutive. We notice that, in a symbolic 
way, w A duj is now a scalar c(x) providing the zero order equation £ r d r c(x) — and the condition 
is c(x) = c = est. The integrability condition (IC) is the Vessiot structure equation: 

wi(<9 2 ^3 - d 3 uj 2 ) + uj 2 (d 3 u;i - <9icj 3 ) + lu 3 (8ilo 2 - d 2 u>i) = c 
involving the only structure constant c. 

For us = (1,— x 3 ,0), we get c = 1. If we choose u> = (1,0,0) leading to c = 0, we may define 
f = {/ e aut(X)\j 1 (f)^ 1 (id) — lj} with infinitesimal transformations satisfying the involutive 
system: 

d 3 e + d 2 e - d^ 1 = o, d 3 e = o, o 2 e = o 

with again 2 equations of class 3 and 1 equation of class 2. The equivalence problem ji(f)^ 1 (uj) = uj 
cannot be solved even locally because this system cannot have any invertible solution. Indeed, 
studying the system ji^g) -1 ^) = co, we have to solve: 

dq 1 ndq 1 „ dg 1 „ dq 1 „ dg 1 „ dg 1 

— — h y = 0, = => — — = 0, — — = 0, — — = 

dy 2 dy 1 ' dy 3 dy 1 ' dy 2 ' dy 3 

by using crossed derivatives. 

EXAMPLE 1.4: (Unimodular contact structure) With similar notations, let us again set a = 
dx 1 ~x 3 dx 2 => da = dx 2 A dec 3 but let us now consider the new Lie pseudogroup of transformations 
preserving a and thus da too, that is preserving the mixed object to — (a, (3) made up by a 1-form 
a and a 2-form (3 with 7 = a A/3 7^ and da = f3 =>• d/3 = 0. Then T is a Lie subpseudogroup of the 
one just considered in the previous example and the corresponding infinitesimal transformations 
now satisfy the involutive system: 

d^ 1 = 0, d^ 2 = 0, d^ 3 = 0, d 2 e + x 3 d 3 e - e = o, o 2 e + o 3 e = o, o 3 e - x 3 & 3 e = o 

with 3 equations of class 3, 2 equations of class 2 and 1 equation of class 1 if we exchange x 1 with 
x 3 , a result leading now to 4 CC. 

More generally, when oj — (a, (3) where a is a 1-form and (3 is a 2-form satifying a A /3 7^ 0, we 
may stedy i/ie same problem as before for the general system = 0, £(£)/? = where £(^) is 

the standard Lie derivative of forms with respect to the vector field £, that is £(£) = i(^)d + di(£) 
if i(^) is the interior multiplication of a form by the vector field £. We now provide details on 



4 



the tedious computation involved as it is at this point that computer algebra may be used. With 
a = (ai),/3 = {fiij = —fiji) we may suppose, with no loss of generality, that ai ^ 0,p23 ^ in 
such a way that aA/3^0<^ aip23 + a 2 (3 3 i + a 3 (3i 2 ^ 0. We can then solve the general equations 
as in the special situation already considered with respect to the 6 leading principal derivatives 
pri = {di^ 1 , <9i£ 2 , <9i£ 3 , d 2 £ 1 , d 2 ^ 2 , d^ 1 } as involution is an intrinsic local property and obtain: 

d^ 1 + + -9ie 3 + ... = o 

Oil Oi\ 

o^e + ^(die + d 3 e) + + ... = o 

P32 P32 

P32 P32 

^l + ^2 + ^3 + =Q 
Oil Oil 

P23 P23 

o 3 e + ^a 3 e + ^d 3 e + ... = o 

Oil Oil 

Solving wth respect to the 6 principal derivatives in a triangular way while introducing the 3 para- 
metric derivatives par = {d 2 ^ 3 7 d 3 £, 2 , d 3 £ 3 }, we obtain for example: 

(i - ^)d 2 e - ^d 3 e + (i - ^)d 3 e - ^o 2 e + ... = o 

P23 Oil P23 Oil P23 Oil p 23 Oil 

Setting now /3 = da and identifying the corresponding coefficients, then £(£)j3 is a linear combi- 
nation of £(£)a and C{£)(5 if and only if we have for example: 

013 03 013 «3 

023 Ctl _ 023 a l 



I _ PlS C2 -i _ 013 Q-2 

023 ai 1 p 23 ai 



and thus: 



Pl3 _ $13 $23 _ $13 _ $12 
$23 $23 $23 $23 $12 

a result, not evident at first sight, showing that the 2-form da must be proportional to the 2-form 
/J, that is da = c'(x)/3 and thus a Ada — c'{x)a A f3. As a A /3 ^ 0, we must have c'(x) = d = est 
and thus da — c'(3. Similarly, we get d(3 — c"a A (3 and obtain finally the 4 Vessiot structure 
equations da — c 1 (3, dj3 = c"a A f3 involving 2 structure constants c = (c', c"). Contrary to the pre- 
vious situation (but like in the Riemann case !) we notice that we have now 2 structure equations 
not containing any constant (called first kind by Vessiot) and 2 structure equations with the same 
number of different constants (called second kind by Vessiot), namely aAda = c'aA/3, df3 — c"aA/3. 
Finally, closing this system by taking once more the exterior derivative, we get = d 2 a = c'dj3 = 
c'c"aAf3 and thus the unexpected purely algebraic Jacobi condition e'e" — 0. For the special choice 
oj = (dx 1 — x 3 dx 2 , dx 2 A dx 3 ) we get c = (1, 0), for the second special choice oi) = (dx 1 , dx 2 A dx 3 ) 
we get c = (0,0) and for the third special choice u) = ((l/x 1 )dx 1 ,x 1 dx 2 A dx 3 ) we get c = (0, 1) 
with similar comments as before for the possibility to solve the corresponding equivalence problems. 



REMARK 1.5: Comparing the various Vessiot structure equations containing structure constants 
that we have just presented and that we recall below in a symbolic way, we notice that these 
structure constants are absolutely on equal footing though they have nothing to do with any Lie 
algebra. 

duj — du = c u u 
$7 — $7 + 77 — 77 = c (5oj — 5uj) 
to A (dco — du) = c 

{da — da — c' (3 
Op = c" aA/3 
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Accordingly, the fact that the ones appearing in the MC equations are related to a Lie algebra is 
a pure coincidence and we may even say that the Cartan structure equations have nothing to do 
with the Vessiot structure equations. Also, as their factors are either constant, linear or quadratic, 
any identification of the quadratic terms appearing in the Riemann tensor with the quadratic terms 
appearing in the MC equations is definitively not correct even though most of mathematical physics 
today is based on such a confusion ([43]). Meanwhile, we understand why the torsion is automati- 
cally combined with curvature in the Cartan structure equations but totally absent from the Vessiot 
structure equations, even though the underlying group (translations + rotations) is the same. In 
addition, despite the prophetic comments of the italian mathematician Ugo Amaldi in 1909 ([1]), 
we do believe that it has been a pity that Cartan deliberately ignored the work of Vessiot at the 
beginning of the last century and that the things did not improve afterwards in the eighties with 
Spencer and coworkers (Compare MR 720863 (85m:12004) and MR 954613 (90e:58166)). 

In the second section of this paper, which is an extended version of a lecture given at the 
international conference SCA (Symbolic Computation and its Applications) 2012 held in Aachen 
(Aix-la-Chapelle), Germany, May 17-20, 2012, we shall recall the definition of the Chevalley- 
Eilenberg cohomology existing for Lie algebras and describe its use in the study of the deformation 
theory of Lie algebras ([28], [29], [45], [49], [50]). We insist on the fact that the challenge solved in the 
remaining of the paper is not to generalize this result to arbitrary Lie equations but rather to work 
out a general framework that will provide exactly the deformation cohomology of Lie algebras in 
the particular case of Example 1. 

In the third section we study the " Vessiot structure equations" in the nonlinear framework and 
explain why they must contain "structure constants" satifying algebraic Jacobi-like conditions. 

In the fourth section, we present the minimum amount of differential geometry (jet theory, 
Spencer operator, (5-cohomology, differential sequences) needed in order to achieve the formal con- 
structions done in the next sections. 

In the fifth and longest section we present for the first time a rather self-contained but com- 
plete general procedure for exhibiting a deformation cohomology for any system of transitive Lie 
equations, even of infinite type as in the case of the (unimodular) contact structure. A link with 
the difficult concept of a truncated Lie algebra will also be provided (Compare to [14], II). 

Finally, we conclude the paper with a few explicit computations based on the previous examples 
while showing out the possibility to use computer algebra techniques. 



2 DEFORMATION THEORY OF LIE ALGEBRAS: 

Let V be finite dimensional vector space over a field k containing the field Q as a subficld and 
set V* = homk{V, k). We shall denote the elements of V by X, Y, Z, ... with components X = (X p ) 
for p = 1, dim(V). 

DEFINITION 2.1: A Lie algebra Q = (V,c) is an algebraic structure on V defined over k by a 
bilinear map [ ] : A 2 V — > V called bracket through the formula ([X, Y]) T = <? X p Y a where the 
structure constants c G A 2 V* ® V are satisfying the Jacobi relations: 

a result leading to the Jacobi identity for the bracket, namely: 

[A, [V, Z\\ + [Y, [Z, X}} + [Z, [X, Y]]=0 , VA, V, Z G V 

Considering a Lie algebra as a point c of the algebraic set defined by the quadratic equations 
J(c) = 0, we may state: 
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DEFINITION 2.2: A deformation ct of c is a curve passing through c in this algebraic set, that 
is to say a set of points ct indexed by a parameter t and such that J(ct) = with ct = c + tC + .... 
As a byproduct, the Lie algebra Qt = (V, c t ) is called a deformation of Q — (V, c). 

Of course, a central problem in the theory of Lie algebras is to exhibit properties of Q that do 
not depend on the basis chosen for V. In particular, if a e aut(V) C V* <S> V, we may define an 
equivalence relation among the structure constants as follows: 

DEFINITION 2.3: c - c ac(X,Y) = c(aX,aY),VX,Y e V or cquivalently c = aaa^c in a 
symbolic tensorial notation showing that c transforms like a (2-covariant, l-contravariant)-tensor 
under the action of a. In particular, we get at once c ~ c <=> J(c) = whenever J(c) = 0. 

When at € aw£(V) is such that ao = idy, then the set of points c* ~ c is a deformation of c 
called trivial and we may state: 

DEFINITION 2.4: a e aut(V) is called an automorphism of <? if c = c and A e end(U) = 
V* (8> V is called a derivation of if A[A,F] = [AY,Y] + [X,AY],MX,Y e V. In particular, if 
a t — a + tA + ... e aut(Q) is such that ao = idy , then A = ^ | t=0 is a derivation of C7. Moreover, 
if we define the adjoint action of V on V by ad(A)Y = [A, Y], then it follows from the Jacobi 
identity that ad(X) is a derivation of called inner derivation. 

DEFINITION 2.5: A Lie algebra Q is said to be rigid if it cannot admit a deformation which is 
not trivial. 

As the algebraic set of structure constants may have very bad local properties, it may be inter- 
esting to study infinitesimal deformations when t is a small parameter, that is t <C 1. In particular, 
we should have dJ g ^ C = for any deformation while a trivial deformation should lead to: 

{cy pa = {AY p (c)i a + {Ay a {cy plx - (A)i(cy pa 

with no way to unify all these technical formulas. 

By chance, such an infinitesimal study can be made easy by means of the Chevalley-Eilenberg 
cohomology that we now describe.. For this, by analogy with the exterior derivative, let us define 
an application d : A r V* ® V — > A r+1 ^* ® V depending on c by the formula: 

d/(X 1; X r+l ) = E^ji-iy+jfdXi, Xj],Xi, Xi, Xj, X r+l ) 

+ Ei(-^ i+1 [Xi,f(x 1 ,...,x i ,...,x r+1 )} 

where a "hat" indicates an omission. Using the Jacobi identity for the bracket and a straightfor- 
ward but tedious computation left to the reader as an exercise, one can prove that dod — 0. Hence 
we may define in the usual way coboundaries B r {Q) as images of d, cocycles Z r (Q) as kernels of d 
and cohomology groups H r (Q) = Z r (Q) / B r (Q) in such a way that B r (Q) C Z r (Q) C A r V* <g> V. 

LEMMA 2.6: C = ^| t=0 S Z 2 {G) and c t - c ^ C = ^| t=0 £ B 2 (^) . Accordingly, a sufficient 
condition of rigidity is H 2 (Q) = 0. 

LEMMA 2.7: A is a derivation ^ A E Z X (Q) and A is an inner derivation ^ A e B^Q). Ac- 
cordingly, the vector space of derivations of Q modulo the vector space of inner derivations is H\{Q). 

Up to the moment we have only been looking at infinitesimal deformations. However, a Lie 
algebra Q may be rigid on the finite level even if H2{Q) ^ 0, that is even if it can be deformed on 
the infinitesimal level ([28]) and we may state: 

DEFINITION 2.8: An element C G Z 2 {Q) is said to be integrable if there exists a deformation 
c t of c such that ^§f | t=0 = ^ ■ ^ i s sa id to be formally integrable if there exists a formal power 
series c t = T7T C » with C o = c, C\ = C. 
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EXAMPLE 2.9: Even if Example 1.4 has nothing to do at first sight with any Lie algebra as we 
already said, we may adapt the previous arguments to c'c" = 0. Indeed, we must have: 

= (c' + tC[ + jC 2 + ...)(c" + tC'l + l ^C'i + ...) - t(c"C[ + c'C'l) + ^(c"C 2 + 2C[C'{ + c'C 2 ') + ... 

As cocycles are defined by the condition c"C + c'C" = 0, it follows that the cocycle (1, 1) cannot 
be integrable at c = (0, 0) as we should have C'C" = 0. 

In order to study the formal integrability of a cocycle C € Z 2 (Q), we shall use a trick by intro- 
ducing two parameters s and t. Then c s+t becomes a deformation of c* in s and we have: 

where one may notice the change C — > X)^Lo Tr^Wi- 

PROPOSITION 2.10: A sufficient condition for the formal integrability of any cocycle is 

#3(0) =0. 

Proof. By definition we have d = d(c) and we may set d(t) = d(c t ) in order to obtain: 

°° t v 

dCi = =► ^ — C„+i = 

L/=0 ^' 

Differentiating with respect to t while using the Leibnitz formula, we get: 
or, in an equivalent way: 

in.. . , j. 

A! M ! V di A 

As the first term on the left belongs to B^Q), it just remains to prove that the sum on the right , 
which only depends on C\, C v , is in Z%{Q). Otherwise, as shown by the previous example, the 
study of these inductively related conditions may sometimes bring obstructions to the deformation 
at a certain order. 

Now we have d o d = => d(t) o d(t) = and thus again: 

A,^i>0,A+/j=j/ ^ 

Accordingly, we obtain successively by using this identity (care to the minus sign): 

\>0,V>0,\+fJ, = l> ^ \>0,fJ,>0,\ + fJ, = L> ^ 

A>0,m>0,A+^=i> ^ a>0,/3>0,a+/3=\ ' 

_^ v\ d°d(t) v (i/-a)l d^d(t) 

a=l y ' ^>0,/3>0,/3+p=i/-a ' ^ 

because we may take into account the integrability conditions up to order v — 1 as a > 0. 

Q.E.D. 

Studying the successive integrability conditions, we get: 
v = dd = 
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and we have the following result not evident at first sight: 

COROLLARY 2.11: The Hessian of the Jacobi conditions provides a quadratic map: 

d 2 1(c) 

H 2 (Q) H Z {Q) : C -> -^-CC 
ococ 

Proof. The previous proposition proves that this map takes C £ ^(5) to Z%(Q) and we just need 
to prove that it also takes C G B 2 (G) to B 3 (Q). For this, let A e V* (g)V be such that C = dA. 
We have: 

-(-^)uc - -Hr)i« ° dA - d ° (-sr)u^ 

as we wished. 

Q.E.D. 



3 VESSIOT STRUCTURE EQUATIONS : 

If X is a manifold, we denote as usual by T = T(X) the tangent bundle of X, by T* = T*(X) 
the cotangent bundle, by A r T* the bundle of r-forms and by 5 g T* the bundle of q- symmetric ten- 
sors. More generally, let £ be a fibered manifold, that is a manifold with local coordinates (x l ,y k ) 
for i = 1, n and k = 1, m simply denoted by (x, j/), projection w : £ — > X : (or, y) — > (or) and 
changes of local coordinates x = <p(ar), y = il>(x, y). If £ and T are two fibered manifolds over X 
with respective local coordinates (x,y) and (x, z), we denote by fxxJ the fibered product of £ 
and J 7 over X as the new fibered manifold over X with local coordinates (x, y, z). We denote by 
/ : X — > £ : (x) — > (x, y = /(x)) a global section of £, that is a map such that n o / = zc?x but 
local sections over an open set U <Z X may also be considered when needed. Under a change of 
coordinates, a section transforms like f(<p(x)) — ip(x,f(x)) and the derivatives transform like: 

||^(x))^(x) = %{x,f{x)) + /(*))*/*(«) 
We may introduce new coordinates (x l ,y k ,y k ) transforming like: 

yld iip r (x) = ^(x,y) + ^(x,y)y* 

We shall denote by J q {£) the g-jei bundle of £ with local coordinates (x l ,y k ,y k ,y k p ...) = (x,y q ) 
called jet coordinates and sections / 9 : (x) — > (x, f k (x), f k (x), f k j{x), ...) = (x,f q (x)) transforming 
like the sections j q (f) : (x) — > (x,f k (x),dif h (x),dijf k (x),...) = (x,j q (f)(x)) where both / g and 
are over the section / of £. Of course J q {£) is a fibered manifold over X with projection 7r g 
while J g+r (£) is a fibered manifold over J 9 (£) with projection 7r^ +r ,W > 0. 

DEFINITION 3.1: A system of order q on £ is a fibered submanifold lZ q C Jg(£) and a solution 
of 7£ 9 is a section / of £ such that is a section of 7£ 9 . 

DEFINITION 3.2: When the changes of coordinates have the linear form x = cp(x), y — A(x)y, 
we say that £ is a vector bundle over X and denote for simplicity a vector bundle and its 
set of sections by the same capital letter E. When the changes of coordinates have the form 
x = <p(x), y = A(x)y + B(x) we say that £ is an affine bundle over X and we define the associated 
vector bundle E over X by the local coordinates (x, v) changing like x = <p(x), v = A(x)v. 

DEFINITION 3.3: If the tangent bundle T(£) has local coordinates (x,y,u,v) changing like 
vP = di^(x)u l , v l = ^r(x,y)u l + §^pr(x,y)v k , we may introduce the vertical bundle V{£) C T{£) 
as a vector bundle over £ with local coordinates (x, y, v) obtained by setting u — and changes 
v l — |pr(x,y)u fc . Of course, when £ is an affine bundle over X with associated vector bundle E 
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over X, we have V(£ ) = £ x x E. 

For a later use, if £ is a fibered manifold over X and / is a section of £, we denote by 
the reciprocal image of V(£) by / as the vector bundle over X obtained when replacing (x, y, v) by 
(x, f(x), v) in each chart. A similar construction may also be done for any affine bundle over £ . 

We now recall a few basic geometric concepts that will be constantly used through this pa- 
per. First of all, if ^, 77 G T, we define their bracket [£,77] E T by the local formula ([£, n]y(x) = 
C(x)d r rj' l (x)-r] s (x)d s C(x) leading to the Jacobi identity [£, [77, (]] + [n, [C,£]] + [C> [?,??]] = 0,V£,?7,C E 
T allowing to define a Lie algebra and to the useful formula [T(/)(£), T(f)(rj)] = T(/)([£, 77]) where 
T(f) : T(X) T(Y) is the tangent mapping of a map / : X ->• Y. 

When I = {ii < ... < i r } is a multi-index, we may set dx 1 = dx 11 A ... A dx lr for describing 
A r T* and introduce the exterior derivative d : A r T* — > A r+1 T* : lu = coidx 1 duj = diUJidx 1 Adx 1 
with d 2 = d o d = in the Poincare sequence: 

/\0rp* d y^iy* j ^ A 2 T* - ) - > A n T* >■ 

The Lie derivative of an r-form with respect to a vector field £ £ T is the linear first order 
operator £(£) linearly depending on ji(£) and uniquely defined by the following three properties: 

1) C(Of = U = Cdif,Vf e A°T* = C°°(X). 

2) C(Z)d - dC(0- 

3) £(f)(a A/9) = A/3 + a A (£(0/3),Va,/3 e AT*. 

It can be proved that £(£) = + where is the interior multiplication (i(£ > )u>)i 1 j r = 
C-'u:...<, and that [£(0, £(0 o £(r?) - C( V ) o £(0 = £([^, r?]),VC, »y G T. 

We now turn to group theory and start with two basic definitions: 

Let G be a Lie group, that is another manifold with local coordinates a = {a 1 , ...,a p ) called 
parameters, a composition G x G — > G : (a,b) — !> ab, an inverse G — > G : a — > a -1 and an identity 
e E G satisfying: 

(a6)c = a (be) = abc, aa" 1 — a~ 1 a = e, ae = ea = a, Va, b,c E G 

DEFINITION 3.4: G is said to act on X if there is a map X x G ->• X : (a;, a) ->• y = aa; = f(x, a) 
such that (ab)a; = a(bx) = abx,\/a,b E G,Vx E X and we shall say that we have a Lie group of 
transformations of X. In order to simplify the notations, we shall use global notations even if only 
local actions are existing. The set G x — {a E G | ax — x} is called the isotropy subgroup of G at 
x E X and the action is said to be effective if ax = x, Va; E X => a = e. 

DEFINITION 3.5: A Lie pseudogroup of transformations T C aut(X) is a group of transfor- 
mations solutions of a system of OD or PD equations such that, if y — f{x) and z = g(y) are 
two solutions, called finite transformations, that can be composed, then z = g o f(x) = h{x) and 
x = ,f^ 1 (y) = g(y) are also solutions while y — x is a solution and we shall set id q = j q (id). 

It becomes clear from Examples 1.1 and 1.2 that Lie groups of transformations arc particular 
cases of Lie pseudogroups of transformations as the system defining the finite transformations can 
be obtained by eliminating the parameters among the equations y q = j q (f)(x,a) when q is large 
enough. The underlying system may be nonlinear and of high order as we have seen. We shall 
speak of an algebraic pseudogroup when the system is defined by differential polynomials that is 
polynomials in the derivatives. Looking for transformations "close" to the identity, that is setting 
y = x + t£(x) + ... when t 1 is a small constant parameter and passing to the limit t — > 0, we 
may linearize the above nonlinear system of finite Lie equations in order to obtain a linear system 
of infinitesimal Lie equations of the same order for vector fields. Such a system has the property 
that, if £, 77 are two solutions, then [£,77] is also a solution. Accordingly, the set 8 C T of solu- 
tions of this new system satisfies [0,0] C and can therefore be considered as the Lie algebra of T. 
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We now turn to the theory proposed by Vessiot in 1903 ([51]) and sketch in a few succes- 
sive steps the main results that we have obtained in many books ([36], [37], [38], [39]). We invite the 
reader to follow the procedure on each of the examples provided for this purpose in the introduction. 

1) If £ = X x X, we shall denote by U q — H q (X,X) the open subfibercd manifold of 
J q (X x X) defined independently of the coordinate system by det(y k ) ^ with source projec- 
tion a q : U q — >• X : (x,y q ) — > (x) and target projection (i q : U q — > X : (x 7 y q ) — > (y). We shall 
sometimes introduce a copy Y of X with local coordinates (y) in order to avoid any confusion 
between the source and the target manifolds. Let us start with a Lie pseudogroup L C aut(X) 
defined by a system lZ q C U q of order q. In all the sequel we shall suppose that the system is 
involutive (see next section) and that L is transitive that is Vx, y e X, 3f e L, y = f(x) or, equiv- 
alently, the map (a q ,(3 q ) : lZ q — > X x X : (x, y q ) — >• (x,y) is surjective. 

2) The Lie algebra G C T of infinitesimal transformations is then obtained by linearization, 
setting y = x + t£(x) + ... and passing to the limit t — » in order to obtain the linear involutive 
system R q = id^iyiJZq)) C J q {T) by reciprocal image with = {£ e T\j g (£) e R q } while taking 
into account the fact that T = id^ 1 (V(X x X)). From now on we shall suppose that R q is transi- 
tive, that is to say the canonical projection 7Tq : J q (T) — > T induces an epimorphism 7Tq : R q — > T 

with kernel R q C R q and we have the useful short exact sequence — > R q — > R q T — ► 0. 

3) Passing from source to target, we may prolong the vertical infinitesimal transformations 
V = V k (y)-^k to the jet coordinates up to order q in order to obtain: 

„*(«)— + + ( ^ y r y s + tot r. J_ 

' yU, dy k dy rUl dy\ ^dy^dy^^ 3 dy^ 10 ' dy% 

where we have replaced j q (f)(x) by y q , each component beeing the "formal" derivative of the 
previous one obtained by introducing di — di + y^+^-^F- 

4) As [9,6] C 9, we may use the Frobcnius theorem in order to find a generating fundamental 
set of differential invariants {& T (y q )} up to order q which are such that $ T (y g ) = $ T (y q ) by using 
the chain rule for derivatives whenever y = g(y) e T acting now on f. Of course, in actual practice 
one must use sections of R q instead of solutions but it is only in section 5 through Definition 5.2 
that we shall see why the use of the Spencer operator will be crucial for this purpose. Specializing 
the $ r at id q (x) we obtain the Lie form & T {y q ) = uj t (x) of lZ q . Finally, if $ T is any differential 
invariant at the order q, then di$ T is a differential invariant at order q+ l,Vi = 1, ...,n. 

5) The main discovery of Vessiot, fifty years in advance, has been to notice that the prolongation 
at order q of any horizontal vector field £ = (a;) commutes with the prolongation at order q 
of any vertical vector field 77 = ^(y)-^, exchanging therefore the differential invariants. Keeping 
in mind the well known property of the Jacobian determinant while passing to the finite point of 
view, any (local) transformation y = f{x) can be lifted to a (local) transformation of the differen- 
tial invariants between themselves of the form u — > \(u,j q (f)(x)) allowing to introduce a natural 
bundle T over X by patching changes of coordinates x = <p(x),u — \(u,j q (ip)(x)). A section lo of 
T is called a geometric object or structure on X and transforms like u>(f(x)) — X(co(x),j q (f)(x)) or 
simply to — j q (f)(uo). This is a way to generalize vectors and tensors (q = 1) or even connections 
(q = 2). As a byproduct we have r = {/ e aut(X)\$ UJ (j q (f)) = j q (f)~ 1 (u) = lo} as a new way 
to write out the Lie form and we may say that V preserves lo. Replacing j q {f) by f q , we also 
obtain lZ q = {f q £ II g |/~ 1 (a;) = lo}. Coming back to the infinitesimal point of view and setting 
ft = exp(t£) G aut(X), \/£ e T, we may define the ordinary Lie derivative with value in a; _1 (F(J 7 )) 
by the formula : 

vi = iu = £(0<" = ^(/*) _1 Mlt=o => = T \ £ (0" = 0} 

We have x — > x — x + t£(x) + ... => u T — >• u T = u T + td^ k Lj,^(u) + ... where /x = (fi\, /x„) is a 
multi-index and we may write down the system of infinitesimal Lie equations in the Medolaghi form: 
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W ee (£(0^Y = -L T k »{u(x))d li S k + Cd r u T (x) = 

as a way to state the invariance of the section uj of J 7 , that is u T — uj t (x) = u T — uj t (x) = 0. 
Finally, replacing j q {£) by a section £ q £ J q {T) over £ e T, we may define i? g C J g (T) ok sections 
by the linear (non-differential) equations: 

0- ee (L(£,)w) T ee -L^(uj(x))^ + C'd r ^{x) = 
and obtain the first prolongation R q +i C J q+ \(T) by adding: 

QT ee (L(^ +1 )ji(u;))[ ee -^(c,(x))^ +li - 9L ^ (X)) 9 ^^ x )^+^ T ^)g+r^(g^ r (^)) = 

6) By analogy with "special" and "general" relativity, we shall call the given section special 
and any other arbitrary section general. The problem is now to study the formal properties of the 
linear system just obtained with coefficients only depending on ji(u>), exactly like we did in the 
examples of the introduction. In particular, if any expression involving uj and its derivatives is a 
scalar object, it must reduce to a constant because T is assumed to be transitive and thus cannot 
be defined by any zero order equation. Now one can prove that the CC for uj, thus for uj too, 
only depend on the $ and take the quasi-linear symbolic form v ee I{u\) ee A(u)u x + B(u) = 
with Mi = (u, u x ), allowing to define an affinc subfibered manifold B\ C Ji{F) over T. Indeed, 
if ^2 A{y q )d x § is a minimum sum of formal derivatives of differential invariants of order q not 
containing any jet coordinate of strict order q + 1, we may suppose by division that the first A 
in the sum is equal to 1. Applying the prolonged distribution of vector fields introduced in the 
step 3 at order q + 1, we obtain a new sum with less terms and a contradiction unless all the A 
are again differential invariants at order q and thus functions of the $ because of the Frobcnius 
theorem. A similar comment can be done for the B. Now, if one has two sections uj and uj of J 7 , the 
equivalence problem is to look for / e aut(X) such that j q (f)^ 1 (uj) = uj. When the two sections 
satisfy the same CC, the problem is sometimes locally possible (Lie groups of transformations, 
Darboux problem in analytical mechanics,...) but sometimes not ([36], p 333). 

7) Instead of the CC for the equivalence problem, let us look for the integrability conditions 
(IC) for the system of infinitesimal Lie equations and suppose that, for the given section, all the 
equations of order q + r are obtained by differentiating r times only the equations of order q, 
then it was claimed by Vessiot ([51] with no proof, see [36], p 313, [39], p 207-212) that such a 
property is held if and only if there is an equivariant section c : T — > T\ : (x, u) — > (x, u,v = c(u)) 
where T\ = Ji(j r )/Bi is a natural vector bundle over T with local coordinates (x, u, v). Moreover, 
any such equivariant section depends linearly on a finite number of constants c called structure 
constants and the IC for the Vessiot structure equations I(ui) = c(u) are of a polynomial form 
J(c) =0. It is important to notice that the form of the Vessiot structure equations is invariant 
under any change of coordinate system. In actual practice, this study can be divided into two parts 
according to the following commutative and exact diagram: 



-> 
I 

-> T 

II 

-> T — > 

I 


Indeed, chasing in this diagram, we discover that tt^ +1 : R q +i — > R q is an epimorphism if and only 
if 7T* +1 : R q+1 — > Rq is an epimorphism and 7Tq +1 : R q +\ — V T is also an epimorphism. Looking at 
the form of the corresponding Medolaghi equations L(^ q )uj = and L(£ q+ i)ji{uj) = 0, these two 
conditions respectively bring the Vessiot structure equations of first kind I*(u\) = not invoving 
any structure constant and the Vessiot stucture equations of second kind = c involving 









1 


I 


9q+l 


= 9q+l 


I 


I 




> Rq+1 


1 


I 


R° q 


► Rq 
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the same number of different structure constants. Such a study, only depending now on linear 
algebraic techniques, can be achieved by means of computer algebra ([3], [30]). 
Finally, looking at the formal integrability of the system B\ C Ji(F) defined by the equations 
A(u)u x + B(u) = and their first prolongation A{u)u xx + u x u x + dB g^ u x = 0, the symbols 
only depend on A(u) and we may obtain equations of the form a{u)u x u x + b(u)u x = by eliminat- 
ing the jets of order 2. Using local coordinates (x,u,v = A(u)u x + B{u)) for T\, and substituting 
(u, v) in place of (u, u x ), we obtain equations of the form a(u)vv + (3{u)v + j(u) — 0. As we may 
suppose that c = for the special section, we finally get equations of the form a(u)II + (3(u)I = 
and it only remains to set I(u\) = c(u) in order to get polynomial Jacobi conditions of degree 
< 2 which may not depend on u anymore because these equations are invariant in form under any 
change of coordinates. 

REMARK 3.6: When q = 1, a close examination of the Medolaghi equations and their first pro- 
longation shows at once that we can choose v = A(u)u x and we get homogeneous Jacobi conditions 
of degree 2. Such a result explains the difference existing between Examples 1.1, 1.3, 1.4 (q = 1) 
and Example 1.2 (q = 2). 



4 LINEAR AND NONLINEAR JANET SEQUENCES: 

Let fi — (hi, ...,/x„) be a multi-index with length = fi\ + ... + /i„, class i if \i\ = ... = = 
^ and fi + \ i = (fj, 1 , ^ + l,fj, i+1 , ...,/x n ). We set y q = {y k \l < k < m, < \/j,\ < q} 
with y k = y k when \n\ = 0. If E is a vector bundle over X with local coordinates (x l ,y k ) for 
i = l,...,n and fc = l,...,m, we denote by J q (E) the g-jet bundle of E with local coordinates 
simply denoted by (x,y q ) and sections f q : (x) — > (x, f k (x), f k (x), f^(x), ...) transforming like 
the section : (.t) — >■ (x,f k (x),dif k (x),dijf k (x),...) when / is an arbitrary section of E. 

Then both / q G J q (E) and G ^(-E 1 ) ar e over f £ E and the Spencer operator just allows 

to distinguish them by introducing a kind of "difference" through the operator £> : J q+ i(E) — > 
T* ® : -> - with local components (9,/ fc (x) - f k (x),d l f k (x) - f k (x), ...) 

and more generally (D/ g +i)^ |f (a;) = dif k (x) — Z^+i^x). In a symbolic way, w/ie« changes of 
coordinates are not involved, it is sometimes useful to write down the components of D in the 
form di — di — Si and the restriction of D to the kernel S q +iT* (g> E of the canonical projection 
: -/g+iCE) -> ^g(-E) is minus the Spencer map 5 — dx i A 5i : S q+1 T* ® E ->• T* ® 5 g T* ® £. 
The kernel of £> is made by sections such that f q+ i = ji(f q ) = 32{f q -i) = ■■■ = j q +i(f)- Fi- 
nally, if R q C Jq{E) is a system of order q on E locally defined by linear equations <fr T (x,y q ) = 
a 1 k IJ '{x)y k = and local coordinates (x, z) for the parametric jets up to order q, the r- prolongation 
R q +r = Pr(Rq) = J r (Rq) H J q+r (E) C J r (J q (E)) is locally defined when r = 1 by the lin- 
ear equations $ r (x,j/ 9 ) = 0, di$ T (x, y q +i) = a^ M (x)y* +1 . + diO T k ll (x)y k l = and has symbol 
9q+r = Rq+r H S q + r T* (g) E C J q + r (E) if one looks at the iop orrfer terms. If G is 
over f q £ R q , differentiating the identity a^(x)f k (x) = with respect to x l and substracting the 
identity a T k f "(x)f k +u (x)+d,a T k f "(x)f k (x) = 0, we obtain the identity a£" (*)($/* (s) (x)) = 
and thus the restriction Z) : — > T* ® R q . This first order operator induces, up to sign, the 

purely algebraic monomorphism — > g q+ \ Ar*®j, on the symbol level ([36], [40], [49]). 

DEFINITION 4.1: R a is said to be formally integrable when the restriction Tr q + r ~^ : R q -\- r +i — ^ 
Rq+r is an epimorphism Vr > or, equivalently, when all the equations of order q + r are obtained 
by r prolongations only, Vr > 0. In that case, R q +\ C Ji(R q ) is a canonical equivalent formally 
integrable first order system on R q with no zero order equations, called the Spencer form. 

DEFINITION 4.2: R q is said to be involutive when it is formally integrable and the symbol 

g q is involutive, that is all the sequences ... A A S T* ® g q+r A ... are exact V0 < s < n,Vr > 0. 
Equivalently, using a linear change of local coordinates if necessary, we may successively solve the 
maximum number /3™, /3™" 1 , (3 q of equations with respect to the principal jet coordinates of 

strict order q and class n, n — 1, 1 in order to introduce the characters a q — m ( q ^^(n-i)\ ~ &l 
for i = 1, ...,n with a™ = a. Then i? g is involutive if R q +i is obtained by only prolonging the (3 q 
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equations of class i with respect to d\, di for i = 1, n. In that case dim(g q +i) — + ... + a™ 
and one can exhibit the Hilbert polynomial dim(R q + r ) in r with leading term (a/n!)r™ when a 7^ 0. 
Such a prolongation procedure allows to compute in a unique way the principal (pri) jets from the 
parametric (par) other ones ([23], [46]). 

REMARK 4.3 : This definition may also be applied to nonlinear systems as well while using a 
generic linearization by means of vertical bundles. For example, with m = 1, n = 2 and q = 2, 
the nonlinear system y 2 2 — |(yn) 3 = 0, j/12 — |(yn) 2 = is involutive but the nonlinear system 
2/22 - 5(yn) 2 = 0, yi2 - yu = is not involutive. 

When i? 9 is involutive, the linear differential operator V : E J q (E) A J q (E)/R q = F of 
order g with space of solutions 6 C E is said to be involutive and one has the canonical linear 
Janet sequence ([39], p 144): 

where each other operator is first order involutive and generates the compatibility conditions (CC) 
of the preceding one. As the Janet sequence can be cut at any place, the numbering of the Janet 
bundles has nothing to do with that of the Poincare sequence, contrary to what many people believe. 

Equivalently, we have the involutive first Spencer operator D\ : Cq — R q Ji(R q ) — > 
Ji(Rq) / Rq+i ~T*<E) R q /8(g q+ i) = C\ of order one induced by D : R q+1 ->• T* <g> R q . Introducing 
the Spencer bundles C r — AT* <g) R q /5(A r ~ 1 T* <g) g q +i), the first order involutive (r + 1)- Spencer 
operator D r+1 : C' r ->• C' r+1 is induced by D : A r T* <g> R q+1 -> A r+1 T* <g> i? 9 : a ® -> 
da <E> £ 9 + (— l) r a A -D£ q +i and we obtain the canonical linear Spencer sequence ([39], p 150): 

as the Janet sequence for the first order involutive system R q +\ C Ji(R q ). Introducing the other 
Spencer bundles C r (E) = A r T* ® J q (E)/S(A r ~ 1 T* ® S g+ iT* ® £) with C r C C r (£), the linear 
Spencer sequence is induced by the linear hybrid sequence: 

O^E^ C (E) A d(E) ^C 2 ^ ... ^ C n — ► 

which is at the same time the Janet sequence for j g and the Spencer sequence for J q+ \(E) C 
Ji(J 9 (£)) ([39], p 153): 

We have the following commutative and exact diagram allowing to relate the Spencer bundles 
C r and C r (E) to the Janet bundles F r = AT* <8> F /c5(A r ~T* ® fti) if we start with the short 



exact sequence — > g q+ \ — > S q+ iT* ig 


£ -> ft 


1 ->■ where fti C T* ( 


8T : 













I 




1 




A r - 1 T*®g q+1 




AT* <g> i? 9 — ► 


C r 


I 




1 


I 


A r ~T* ® S q+ iT* (g> E 




A r T*®J q (E) — ► 


C r (E) 










A r -T*®/n 


f ) 


AT* ® F — ► 


F r 






1 


I 














In this diagram, only depending on the linear differential operator T> = $ o j ?) the epimorhisms 
$ r : C r (E) — > F r for < r < n arc induced by the canonical projection $ = $ : Cq(E) = 
J q (E) — > J q (E)/R q = Fo if we start with the knowledge of R q C J q (E) or from the knowledge of 
an epimorphism $ : J q (E) — > F if wc set R q = fcer($). It follows that the hybrid sequence projects 
onto the Janet sequence and that the kernel of this projection is the Spencer sequence. Also, chas- 
ing in the diagram, we may finally define the Janet bundles, up to an isomorphism, by the formula: 
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F r = A r T* ® J q (E)/(A r T* ®R q + 6{A r ~ 1 T* ® S 9+1 T* ® Ej) 

that will be crucially used in the next section dealing with the deformation theory of Lie equations 
when E = T , R q C J q (T) is a transitive involutive system of infinitesimal Lie equations of order q 
and the corresponding operator V is a Lie operator. 

DEFINITION 4.4: The Janet sequence is said to be locally exact at F r if any local section of F r 
killed by XV+i is the image by V r of a local section of F r _i over a convenient open subset. It is 
called locally exact if it is locally exact at each F r for < r < n. The Poincare sequence is locally 
exact but counterexamples may exist ([36], p 202). 



5 DEFORMATION THEORY OF LIE EQUATIONS: 

Having in mind the application of computer algebra to the local theory of Lie pseudogroups, 
we want first of all to insist on two points which have never been emphasized up to our knowledge. 

In order to motivate the first point, we sketch it on an example. 

EXAMPLE 5.1: Let a — ai(x)dx l be a 1-form as in Example 1.1 or 1.4. If we look for infinites- 
imal transformations preserving a, we have to cancel the Lie derivative as follows: 

(C(0a)i = Or{x)diC + Cdr^ix) = 

As a byproduct, we have a well defined Lie operator T> : T — > T* : £ — > XJ>£ = £(£,)a such 
that X>£ = 0,T>n = => V[(,rj\ = because of the well known property of the Lie derivative 
[£(£), C{r})] = C{£) o C(rf) — C(rf) o C{£) = £([£, rf\) and such a property can be extended to tensors 
or even any geometric object. Accordingly, it is usual to say that, if we have two solutions of the 
system, then their bracket is again a solution. However such a result is coming from mathematics 
and cannot be recognized by means of computer algebra, contrary to what is sometimes claimed. 
Surprisingly, the underlying reason has to do with formal intcgrability. Indeed, if we study the first 
derivative of the bracket [£,77], it involves in fact the second derivatives of £ and n and sometimes 
things may change a lot. For example, if da ^ 0, as C(£)da = di(£)da = dC{^)a = 0, then the first 
order equations brought by C(^)da = may not be linear combinations of the first order equations 
brought by C(£)a = 0. In particular, if n — 2 and a — x 2 dx 1 , then da — —dx 1 A dx 2 and the new 
first order equation di^ 1 + c^C 2 = 0, which is automatically satisfied by any solution of the system 
R\ C Ji{T) defined by C(£)a = 0, is not a linear combination of the equations defining R\, that is 
Ri is not involutive as it is not even formally integrable because R^ = ir 2 (R 2 ) C R\ with a strict 
inclusion. Hence it becomes a challenge to define a kind of bracket for sections of Ri C J\ (T) and 
not for solutions as usual, that is to say independently of formal inteyr ability. This idea, which is 
a crucial one indeed as it will lead to the concept of Lie algebroid, is to replace the classical Lie 
derivative £(£) for any £ € T by a formal Lie derivative L(£i) for any £1 e Ji(T) over ( e Tin such 
a way that £(£) = and to compute the bracket [£(£1), L(r)i)] = L(£i)oL(r)i)—L(r]i)oL(£i) 

in the operator sense in order to be sure that the new bracket on Ji(T) will satisfy the desired 
Jacobi identity. We obtain successively: 

(£(£i)a)i = a r {x)d + C'd r oti{x) = 0, (L(r?i)a)j = a s {x)r) s j + T] s d s aj(x) = 

L([^, m ])a = a k (x)(Cd rV * + - nti k s - n s d s tf) + ([£, n]) k d k a t {x) = 
as a way to define: 

The induced property C Ri can therefore be checked linearly on sections and no longer on 

solutions. In particular, we may exhibit the section j^ 1 = 0,£ 2 = 0,(1 = O,^ = 0, £ 2 = 0,£ 2 = 1} 
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of Rt, even if £ + g ^ 0. 



Now, using the algebraic bracket {j q +i(£), j q +i(i])} = jgd^v])^^ 7 ] G T, we may obtain by 
bilinearity a differential bracket on J g (T) extending the bracket on T: 

which does not depend on the respective lifts and 77^+1 of £ 9 and rj q in J 9+ i(T). Applying j q 
to the Jacobi identity for the ordinary bracket, we obtain: 

{€q+l,{Vq+2,Cq+2}} + {Vq+l,{Cq+2,€q+2}} + {tq+l,{€q+2,Vq+2}} = , Vq+2 , (g+2 £ Jq+2(T) 

As we shall see later on, this bracket on sections satisfies the Jacobi identity and the following 
definition is the only one that can be tested by means of computer algebra: 

DEFINITION 5.2: We say that a vector subbundle R q C J q {T) is a system of infinitesimal 
Lie equations or a Lie algebroid if [R q ,R q ] C R qi that is to say [£ g ,r7 q ] £ R q ,Vt; q ,r] q £ R q . The 
kernel i?° of the projection 7Tq : R q — > T is the isotropy Lie algebra bundle of 1Z® = id~ (1Z q ) and 
[Rq, Rq} C i?° does not contain derivatives, being thus defined fiber by fiber. 

Of course, another difficulty to overcome in this new setting, is that we have no longer an 
identity like dC{£)a — C{£)da = but it is easy to check in local coordinates that: 

(dL(^)a - L(£i)da) tf = a r (x)(d t % - d&) + (d^ r - ^)d r a,(x) - (d 3 C - ^)d r a t (x) 
and the Spencer operator allows to factorize the formula if we notice that: 

{di€ r j -d j G) = (d i c j -G j )-{d j s-G j ) 

We finally obtain: 

*(C(i))i(C(2))(d£(£i)a - m)da) = i(( {2) )L(i(( {1) )DZ 2 )a - i{^)L{i{^ 2) D^)a 
and more generally: 

LEMMA 5.3: When a £ A r_1 T* we have the formula: 

r 

*(c ( i))...i(C(r))(^Ki)« - m)da) = ^(-i) s + i *(c (1 ))... J (C( S ))...«(c w )i(«(c (S )^6)a 

s=l 

which does not depend on the lift £2 G J 2 {T) of £1 € J\{T). 

In order to understand the second point, we have to revisit the work of Vessiot. Indeed, if we 
have a geometric object, that is a section w of a natural bundle T of order q, then we may consider 
the system TZ q = {f q £ U q | Z" 1 ^) = uj} of finite Lie equations and the corresponding linearized 
system R q = {£ q £ J q (T) \ L{^ q )uj = 0} of infinitesimal Lie equations, both with the particular way 
to write them out, namely the Lie form and the Medolaghi form as it becomes clear from Example 
1.1 to 1.4. As a byproduct, when constructing the Janet sequence, we can write F = J q (T)/R q 
but we can also use the isomorphic definition Fq = w (V(.F)) depending on whether we want to 
pay attention to the system or to the object. The main idea of deformation theory will be to begin 
with the second point of view and finish with the first. Starting with a system R q C J q {T), we 
shall suppose that R q is transitive with a short exact sequence —> Rq ^> R q — > T ^ and, 
whatever is the definition of Fq, introduce an epimorphism $ : J q (T) — > Fq while considering the 
following commutative and exact diagram: 
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> 70/71-1 
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II 




— > ^(T) 


— > Fa 


1 
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T 


T 


— > 

















0^ R° ^ J°(T) — ► F — ►() 

I I II 

0^ R a — ► J D (T) A F ^0 

- 



The next definition will also be crucial for our purpose and generalizes the standard definition: 

DEFINITION 5.4: We say that a vector bundle F is associated with i? g if there exists a first 
order differential operator L(£ q ) : F — s> F called formal Lie derivative and such that: 

1) HZq + ri g ) = L(C q ) + L( Vq ) V£ 9 , T] q G R q . 

2) £(/£,) = fL(£ q ) V£ 9 G i? 9 ,V/ G C°°{X). 

3) [L(^),L(7 7 ,)] = L(^)oL(t 7 ,)-L(» 7 ,)oL(^) = L([^,»7,]) V£ 9 ,r/ 9 Gi? 9 . 

4) m q ){fr,) = fL(^)n + (U)V VC 9 G R q ,Vf G C°°{X),Vr, G F where 6/ = i(04f- 

As a byproduct, if E and F are associated with R q , we may set on E ® F: 

Lfe)(r? ® C) = ® C + »7 ® £(£,)C VC 9 G R q , V77 G F, V( G F 

REMARK 5.5: If 6 C T denotes the solutions of R q , then £(£) = L(j q (£,)) is simply called the 
classical Lie derivative but cannot be used in actual practice as we already said because 8 may be 
infinite dimensional as in Example 1.3 and 1.4. We obtain at once: 
l)C(Z + r))=C(0+C(r,) V(,r,ee. 

3) [C(£),C(ri)]=C([ZM) V^GG. 

4) = fC(0v + (U)v ve g e,v/ g c°°(x),v v g f. 

The extension to tensor products is well known. 

The following technical proposition and its corollary will be of constant use later on: 
PROPOSITION 5.6: We have: 

i(()D{Z q+1 ,r] q+1 } = {i(()DZ q+1 ,ri q } + {Z q ,i(()D Vq+1 } 

Proof. We have: 

X-\-LL — l/ 

Now, caring only about we get: 

and the Proposition follows by bilinearity. 

Q.E.D. 

The proof of the following proposition is similar and left to the reader as an exercise: 
PROPOSITION 5.7: We have the formula: 

i(0DK,+i,»7 g +i] = [i{QDi q+u r, q ] + ^ q ,t(0D Vq+1 ]+i(L( m )0D^ q+1 -z(L(^)0D Vq+1 
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COROLLARY 5.8: If R q C J q {T) is such that [R g ,R q ] C R q , then R q+1 C J,+i(T) satisfies 
C Rq+i even if i? g is not formally integrable. 

EXAMPLE 5.9: T and T* both with any tensor bundle are associated with Ji(T). The case 
of T* has been treated at the beginning of this section while for T we may define L{^\)r} = 
[£, rj\ + i{ri)D^ = {a,ii(?7)}- We have indeed ^drV k ~ V s d s £ k + r? s (d s £ fe - & = -rf£ + £ r drV k 
and the four properties of the formal Lie derivative can be checked directly as we did for T* . Of 
course, we find back £(0?7 = [£, 77], V£, 77 £ T. 

More generally, we have in a coherent way: 

PROPOSITION 5.10: J q (T) is associated with J g+ i(T) if we define: 

and thus i? 9 is associated with R q+ \. 

Proof. It is easy to check the properties 1, 2, 4 and it only remains to prove property 3 as follows. 
[L^ g+1 ),L(r] g+1 )]C q = L(^ +1 )({»7, +1 ,C g +i} + *(»7)I>C,+i) -% + i)(fe + i,C,+i} + 

= { ? ?9+2j C9+2}} - {£g+2, Cg+2}} 

+{£ g +i,«fa)-DC 9 +2} - {r) g +i,i{£)DC q+2 } 

+i(0 D iVq+2, Cq+2} - 1 {v)D{£ q+2 , (q+2 } 

+z(0 J D(*(r?) J DC 9 +2) - ifa)^)^^) 

= {U<J+2, 17^+2}, Cq+l} + {i(0Dr]q+2,(q+l} ~ {i(v)D£q+2 , C<Z+1 } 

+z(K,77])£»C 9 +i 
= {[£ 9 +i,»?9+iLCh-i} + *([£, 

by using successively the Jacobi identity for the algebraic bracket and the last proposition. 

Q.E.D. 

COROLLARY 5.11: The differential bracket satisfies the Jacobi identity : 

[tq, [Vq,Cq]] + [Vq, [Cq,tq]] + [C„ [Cq,Vq]] = V£„lfc,C, & Jq{T) 

PROPOSITION 5.12: We have the formula: 

i(C)(DL^ q+2 ) Vq+1 - L(Z q+1 )D Vq+1 ) = L(i(()DZ q+2 )r] q 
Proof. Using Proposition 5.6, we have: 

i{QDL^ q+2 )r, q+1 = i(C)£>{£,+2,»7,+2} + H£)Di(§Dn q+ 2 

= {i(0D^ q+2 ,r, q+1 } + {^ q+1 ,i(C)Dr )q+2 } + i(0Di(0Dr 1q+2 

and we must substract: 

i(0£(£,+i)I>>7,+i = L(t q+1 )(i(()Dr ] q +1 )-i(L(C 1 )C)DT 1 q +1 

= {Z q+u i(0Dr} q+2 } + i(0Di(C)Dr} q+2 -i(L(Z 1 )()Dr )q+1 

in order to obtain for the difference: 

L(t(C)DCq +2 )vq - i{i{0D^)D nq+1 + i(L(^)0Dt?,+i + 7(0^(0^+2 - mDi(()Dr] q+2 
Finally, the last four terms vanish because £(£i)C — *(0-^£i = [£j C] an d: 

i(0Di^)D Vq+ 2-i(0Di(0Dr, q+ 2 = -*([£, CD^+i- 

Q.E.D. 

Combining this proposition and Lemma 5.3, we obtain: 
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PROPOSITION 5.13: When A r q ~\ G /\ r-1 T* ® J q+ i(T), we have the formula: 

r 

i(t {1) )...i(( {r) )(DL(t q+2 )-L(t q+1 )D^ 

where A^" 1 G a^T* ® J q {T) is the projection of A r q ~\. 

Proof. With a G A r_1 T* and r/ g+ i G J g+ i(T), we obtain successively: 

DL{^ q+2 ){a®r) q+1 ) = D(L(£ 1 )a®ri q+ i+ a® L(i q+2 )ri q+ i) 

= dL(^)a ® Vq + {-ly-^L^a) A D Vq+1 

+da ® Lfe+i)^ + (-l) r_1 a A DLfe+sK+i 

L(^ + i)D(Q!®77g + i) = i(^+i)(da«)?7 9 + (-l) r - 1 aAL>?7 9+ i) 
= L(£i)da ® r\ q + da <g> L(£ 9+ i)?7 9 

+(-ir- 1 L(6)a A Dt7 9+1 + (-l)'- 1 a A L(Z q+1 )Dr, g+1 

and obtain y substraction: 

(DL(t q+2 )-L({ q+1 )D)(a®r, q+1 ) = (dL(^)-L(^)d)a^r ]q ^-iy- 1 aA(DL^ q+ 2)-L^ q+1 )D)r lq+1 
and the proposition follows by skewlinearity. 

Q.E.D. 

PROPOSITION 5.14: We have the formula: 

HQiVq, Cq} = {L(tq+l)Vq, Cq} + {Vq, Htq+lKq} 

which does not depend on the lift G J q+ \(T) of £ q G J q (T). 

Proof. Using the Jacobi identity for the algebraic bracket and proposition 5.6, we obtain: 

L (Cq){VqXq} = {tq , {Vq+1 , Cq+1 } + i(€) D {Vq+l , Cq+1 } 
= {{Cq+1, Vq+l}, Cq} + {Vq, {Cq+1 , Cq+1 } } 

+ {i(ODT) q+1 ,Cq} + {Vq, mDCq+l} 
= { L (Cq+l)Vq,Cq} + {Vq,HCq+l)Cq} 

Q.E.D. 

Finally using Proposition 5.12, we obtain at once: 
COROLLARY 5.15: We have the formula: 

L (Cq+l)[Vq,Cq] = [HCq+1 )Vq , Cq] + Vlq , HCq+1 )Cq] 

+L(i(()D£ q+2 )r) q - L(i(r,)D£ q+2 )( q 
which does not depend on the lift £ q+2 G J q+2 (T) of ^ q+ i G J q+ \(T). 

Before going ahead, let us stop for a moment and wonder how we could proceed for generalizing 
the deformation theory of Lie algebras by using the Vessiot structure equations even though we 
know that the structure constants have nothing to do in general with any Lie algebra. Of course we 
could start similarly from the Jacobi relations but, if we do want to exhibit a kind of cohomology, 
we should be able to define a trivial deformation, that is the analogue of a change of basis of the 
underlying vector space V of the Lie algebra Q in such a natural way that it could induce a change 
of the structure constants which is surely not of a tensorial nature anymore. 

The following "trick", already known to Vessiot in 1903 ( [51], p 445), is still ignored today. For 
this, assuming that the natural bundle T is known, let us consider two sections u> and w giving 
rise respectively to the systems R q and R q of infinitesimal Lie equations: 

Rq W EE -LI»(UJ(X))$ + CdrU r (x) = 
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Rq 



sr = - 



L^{Q{x))e tl +Cd r u> T {x)=Q 



and define the following equivalence relation: 

DEFINITION 5.16: Q ~ lo & R q = R q 

The study of such an equivalence relation is not evident at all and we improve earlier presen- 
tations (compare to [36], p 336). First of all, having in mind what we did for Example 1.4, we 
shall use a solved form of the system obtained by choosing principal jets or, equivalently, choosing 
a square submatrix M = (M(u)) of rank dim(F ) in the matrix L = (L(uj) defining R° q with 
dim(Jg(T)) columns which describe the infinitesimal generators of prolongations of changes of 
coordinates on X acting on the fibers of F and dim(F ) — m rows. Of course, a major problem 
will be to obtain intrinsic results not depending on this choice. The columns of L are thus made 
by vector fields L% = L^ m (m)^f that we can therefore separate into two parts, namely the vectors 
L CT = MJ(jt)g^r for a = 1, ...,dim(Fo) and the vectors L m+r = £° n+r (u)L a for r = 1, ...,dim(R^) 
obtained by introducing the stationary functions £{u), also called Grassmann determinants, while 
describing the matrix M^ 1 L = {idp Q ,£(u)). We are therefore led to look for transformations 
u = g(u) of the fibers of T such that: 



In order to study such a system and to prove that it is defining a Lie pseudogroup of transfor- 
mations, let us notice that the first conditions are equivalent to saying that the transformations 
u = g{u) preserve the vector fields L a and also the vector fields L m+r according to the second 
conditions. It follows that the transformations u = g(u) preserves the vector fields L^, a property 
thus not depending on the choice of the principal jets. In addition, we have: 

PROPOSITION 5.17: The Lie pseudogroup of transformations of the fibers of T that we have 
exhibited is in fact a Lie group of transformations, namely the reciprocal of the lie group of trans- 
formations describing the natural structure of T . 

Proof. The defining system is finite type with a zero first order symbol. If W = W T (u)-^ is an 
infinitesimal transformation, we obtain therefore the Lie operator [W, L%] = 0, VI < \p\ < q, Vfc = 
1, . . . , m. Indeed, if W\ and Wi are two solutions, then [W\, W2] is also a solution because of the 
Jacobi identity for the bracket and there are at most dim(Fo) linearly independent such vector 
fields denoted by W a . It follows that [Wi, W2] — Pi 2 ( u )W a and we deduce from the Jacobi identity 
again that ■ Pi 2 (u) = => Pi 2 (u) = pf 2 = cs t because rk(L^) — dim(F ). Accordingly, the 
W a are the infinitesimal generators of a Lie group of transformations of the fibers of T and the 
effective action does not depend on the coordinate system. 



DEFINITION 5.18: These finite transformations will be called label transformations and will 
be noted u — g(u,a) where the number of parameters a is < dim(F ). 

If R q is formally integrable/involutive, then R q = R q is also formally integrable/involutive and 
thus I(ji(uj)) — c(uo) «4> I(ji(u>)) = c(lj) with eventually different structure constants. 

COROLLARY 5.19: Any finite label transformation u — g(u,a) induces a finite transfor- 
mation c = h(c, a) of the structure constants which is not effective in general and we may set 
lD ~ co => c ~ c. 

It now remains to exhibit a deformation cohomology coherent with the above results. 

DEFINITION 5.20: When F is a vector bundle associated with R qi we may define T = T(F) = 
{n G F\L(£ q )ri = 0,V£ g G R q ] and the sub-vector bundle E = {77 G F\L{£, a q )n = 0,V£° G i?"} C F 
in such a way that T C E C F. 



{M- 1 y T {u)du T = {M- 1 y T {u)du T , 



£° l+r (u) = £°{u) 



Q.E.D. 
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In order to look for T in general, we shall decompose this study into two parts, exactly as we did 

in section 3.7, by using a splitting of the short exact sequence 0— > R® R q T ^ called R q - 
connection, namely a map Xq '■ T — > R q such that 7Tq o x q = idr, in order to have R q ~ R q © \q (T) ■ 
Such a procedure does not depend on the choice of Xq because, if x q is another i? g -connection, 
then (x q — Xq){T) G Rq- An i? 9 -connection may also be considered as a section Xq G T* <g> R q 
over ?c?t G T* ® T and X«j — Xg G T* ® i?° in this case. It follows that we have equivalently 
T = {?] G _E|L(x 9 (C)) 7 ? = 0,V£ G T} and we may define a first order operator V : E ^ T* ® E 
with zero symbol, called covariant derivative, by the formula (V • 7/)(£) = V^?7 = L(xq(£))r]. We 
may now extend V to a first order operator V = AT* ® E -> A r+1 T* ® £ by the formula: 

V(a ® ??) = da <g> 77 + (-l) r a A Vj?, Va £ A r T*, Vr/ G £ 

LEMMA 5.21: With V 2 = V o V, we have: 

(V 2 ry)(e,0 = L([ Xq (0,X q (0] - X g ([CA}))v = => V 2 = 

Proof. We have: 

V 2 (a<g>r/) = V(da®r; + (-l) r aAVf;) 

= d 2 a + (-l) r+1 da A Vry + (-l) r da AVi)+ (-l) r a A V 2 r/ 
= (-l) r aAV 2 r7 

Setting V?7 = dx l Vii], we get: 

V 2 r/ = \7{dx I V l n) = dx l A dx j VjViT] = -\^dx l A dar^ViVj- - Vj-Vi)rj 
and we have just to use the fact that V» = L(x q (di)) with [di, <9j] = 0. Indeed, we have: 

[x 9 (0,x 9 (0] -x,(K,H) - {x 9 +i(0,x 9 +i(C} + i(0^x 9+ i(0 -i(£)D Xq+ i(0 -x 9 (K,Cl) 

The first term in the right member is linear in £ and £ while the sum of the others becomes also 
linear in £ and £ because: 

The proposition follows from the fact that xo — idr and E is i?°-invariant. 

Q.E.D. 

Hence, we obtain by linearity the V- sequence: 

— >T — ► E T* <g> E A 2 T* (8) £ ... A"T* ® £ — > 

which does not depend on the choice of the connection and is made by first order involutive oper- 
ators. It follows that T can be locally described by a linear combination with constant coefficients 
of certain sections of E C F and we may therefore set dim(T) = dim(E) < dim(F). The use of 
computer algebra will essentially be to compute these dimensions by using linear algebra combined 
with homological algebra techniques. 

We now provide a few definitions: 

DEFINITION 5.22: When 9 is given, we may define: 
Centralizer C(9) = {n £ T|[£, n] = 0, V£ £ 9}. 

Center Z(Q) = {n £ 9| [£,77]= 0, V£ £ 9}. 

Normalizer N(Q) = {r? £ 77] C 9,V£ G 9}. 

It is essential to notice that these definitions are not very useful at all in actual practice when 
9 is infinite dimensional. 
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PROPOSITION 5.23: C(9) = T(T). 

Proof. If i?i = 7r?(-Rg) C Ji(T), it follows from Example 5.9 that T(T) = {77 G T\L(^i)rj = 0,V£i € 
that is to say {T(T) = {77 G T|{£i, ji(7?)} = 0,V£i G Ri} if we choose ji (77) as a lift of 77 
in Ji(T). In particular, if £ G 9 and thus ji(£) G i2i, we have {ji (C): Ji C 7 ?)} = [^v] an d thus 
T(T) C (7(0). 

Now, j g _i(K,»7]) - {^(O.J'afa)} and thus C(6) = {77 G ^(77)} - 0,V£ 9 G providing 

by projection {€i,ji(r})} = 0, that is C(9) C T(T) and thus C(9) = T(T). 

Q.E.D. 

It follows that Z(Q) = 9nC(9) =4> Z(Q) = {r, G Tpry = 0,L(£i>7 = 0,V£i G Ri} and Z(9) is 
made by sections of T(T) killed by V. The study of N{&) is much more delicate and we first need 
the next proposition where we notice the importance of involution or at least formal integrability. 

PROPOSITION 5.24: The Lie operator T> : T — > Fq induces a homomorphism of Lie algebras 
V : T(T) — > T(i*o) where the bracket on T(T) is induced by the ordinary bracket on T and the 
bracket on T(F ) is induced by the differential bracket on J q {T). 

Proof: We already know that T inherits a structure of Lie algebra on sections from the ordi- 
nary bracket of vector fields and the situation is similar for J q (T) with the differential bracket. 
Now, if L(£i)j? = and L(£i)C = 0, it follows from Corollary 5.15 that L(£i)[?7,C] = and 
[T(T),T(T)] G T(T). 

Similarly, as ^ = J q (T)/R q = J° q {T)/R q , if Lfe)< = G R° q and L(QC° q ^ G R° q , 

then L(^ q )[r] q , C q ] G R q according to the Jacobi identity for the bracket on J q {T) and the fact that 
[R q ,R q ] C Rg. Also, if L({; q+ i)i] q G R q and L(^ q+ i)( q G R q , then it is less evident to prove that 
L(£q + i)[r) q , ( q ] G R q . For this, using Corollary 5.15, if we set L(£ q+ i)i] q = 9 q G R q , it is sufficient 
to notice that [9 q ,(q] = L(9 q+ i)( q — i(QD6 q+ \ G R q because both terms do belong to R q . Thus, 
in any case, we obtain [T(F ), T(F )] G T(F ). 

Finally, we have j q ([r], Q) — \jq(v)i3q(Q] an d we may take j q (r]) as a representative of Vr) in J q (T). 
We shall prove that, if 77 G T(T), that is if £(£1)77 = 0, V£i G then ,7,(77) G J,(T) is such that 
L(£, q +i)j q (v) = 0- Introducing i? 2 = 7r|(i? g ) C pi(R\) and choosing any £ 2 G -R2 over £1 G we 
have 7To(£(£2)ji(?7)) = £(£i)?7 = =>■ £(£2)ji(f?) G T* ® T. However, using Proposition 5.11 and 
the fact that Dj q {r}) = 0, we obtain: 

i(C)£>(L(6)ii(r?))=i(i(C)^2)»7 = 
because i? 2 C Pi(Ri) => DR 2 C T*<g>Ri and thus L{£,2)ji(ji) = because there is a monomorphism 
(even an isomorphism) — > T* ® T A T* T. Supposing by induction that L{^ q )j q -\{r{) = 0, we 
should obtain in the same way: 

i(QD(L(Z q+1 )j q (ri)) = LiiiQD^jq^in) = 

because R q +i C pi(R q ) => DR q+ i c T* ® R q and thus L(^ q+ i)j q (rf) = because there is a 

monomorphism — ► 5,+iT* ®T \t* ® S q T* ® T and the restriction of £> to a symbol is —5. 

Q.E.D. 

As the reader will discover in the last computational section, the study of the normalizer is 
much more delicate. 

DEFINITION 5.25: The normalizer f = N(T) of F in aut(X) is the biggest Lie pseudogroup 
in which Y is normal, that is (roughly) N(Y) = f = {/ G aut(X)\f o f o /- 1 G T, V/ G T} and we 
write T < 7V(r) G atrf(A). 

Cf course, N(Q) will play the part of a Lie algebra for iV(r) exactly like 9 did for F. However, 
we shall see that N(T) may have many components different from the connected component of the 
identity, for example two in the case of the algebraic Lie pseudogroup of contact transformations 
where N(T)/T is isomorphic to the permutation group of two objects, a result not evident at fist 
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sight. Passing to the jets, we get j q (f 0/0/ x ) = j q (f) o j q (f) 1 o j q (f) = lo O 

jq(.fy 1 (jq(fy 1 ( UJ )) = jqify 1 ^), that is to say = w if we set = w and we 

find back the equivalence relation of Definition 5.16. It follows that r = {/ e aui(X)|j g (/) _1 (o;) = 
g(u>,a),h(c,a) = c} is defined by the system 1Z q +i — {f q +i G n 9+ i|/ g+ i(i?q) = R q } with lineariza- 
tion = {£ q+1 \L(£ g+1 )r} q G i? g ,V?7 9 G that is to say + i(£)Dr] q+1 £ fi, O 
{£ 9+ i,?7 9+ i} G Accordingly, the system of infinitesimal Lie equations defining = iV(©) can 
be obtained by purely algebraic techniques from the system defining 0. In particular, we notice 
that 7Tq +1 : R q +i — > T is an epimorphism because 7Tq +1 : R q +i — > T is an epimorphisme by assump- 
tion and R q +\ C R q +\. We obtain on the symbol level {<7 9 +i, ?7<j+i} C g g and thus <5g g +i C T* ®g q 
leading to C g q+1 = Pl (g q ) and thus g q+1 = g q+1 because R q+1 C C <? 9+1 . Us- 
ing arguments from <5-cohomology, it can be proved that R q +i is involutive when R q is involutive 
([36], p 351, 390). Another proof will be given in Corollary 5.52. 

With more details, using the result of Proposition 5.17, we get the following important local result: 
PROPOSITION 5.26: T = T(F ) = {n r (x) = A a W^(u(x))\A = est} 

Proof. Recalling that Fq — w^ 1 (V (J-)) , we shall first study the natural bundle J-q = V(F) of order 
q. Adopting local coordinates (x, u, v), any infinitesimal change of source x = x + tt;(x) + ... can be 

lifted to J"o with u T = u T + t£*(x)L^(u) + v T = v T + 1 dL § u i U ^ v a + according to the definition 
of a vertical bundle provided by Definition 3.3. The corresponding infinitesimal generators on Tq 
will be: 

(' Wjp HW's + -^Vjf) ' S M < ..ve, e J,(T) 

It follows that a section e : J 7 — > Tq : (a;, u) — >• (x,u,v = e(x, u)) will be equivariant, that is 
v — e(x, u) = => v — e(x, u) = if and only if u = e(u) satisfies: 

Hence, we have e] = => e T (w) = A a WT(u) with A = est. 

As another approach, working directly with Fo, we may consider the invariance of the section 
u — uj(x) = 0, v — fl(x) = and get: 

+?dru T {x) = , - dLTk 2^ x)) ^^K+edr^(x) = 

The first condition brings at once £ q G R q and we obtain therefore the following central local result 
for F : 

To = {n g F | - ^^n^ + ra,.r» T = , ve 9 g izj 

As usual, the study of this system can be cut into two parts. First of all, we have to look for: 

E Q = {ne F \ dLT ^} x)) ^e, = 0,< G R° q } C 

In a symbolic way with pri(R q ) = {£ CT } and par(R q ) = {£ m+r }, we get: 

^lf2^ + d ( £ m+r L *) n gm+r = q 

du du 

whenever £ CT + £^ n+r S, m+r = and thus |^f2 = in agrement with Proposition 5.17. 

Then we have to introduce an i? g -connection, that is a section \q G T* ® i? 9 over G T* ® T 

such that: 



23 



and it just remains to study this last system for £l T (x) = A a {x)W^{uj(x)) as it does not depend 
on the choice of the connection \q- Substituting, we obtain successively: 

B1 T ^ 8W T 

As [L%, W a ] = and the action is effective, we finally obtain diA a = that is A = est. 

Q.E.D. 

The following corollary is a direct consequence of the above proposition and explains many 
classical results as we shall see in the last computational section. 

COROLLARY 5.27: We have the relation: 

9 = {£ G T\V£_ = £(f)w = 0}=> N(&) = = £(f)w = AW(w) G T } 

REMARK 5.28: If we set R q = ir^ 1 (R q +i) , then R q +\ is obtained by a procedure with two 
steps. First of all, we obtain R q = {£ 9 e J q {T)\L{^ q )uj = A(x)W(uj) G E } by eliminating the 
infinitesimal parameters A by means of pure linear algebra. Then, we have to take into account 
that A — est in the above corollary and another reason for understanding that g q +i = g q +i though 
we have only g q C g q in general. Such a situation is well known in physics where the Poincare group 
is of codimension one in its normalizer which is the Weyl group obtained by adding dilatations. 

According to the definition of the Janet bundles at the end of section 4, using the inclusion 
AT* ® S q T*®T C AT* ® J°(T) C AT* ® J q (T), in the case of an involutive system R q C J q {T) 
of infinitesimal transitive Lie equations, we have: 

F r = AT* ® J°(T)/(AT* ® R° q + <S(A r_1 T* ® S q+1 T* ® T)) 
As i?^ is associated with R q and J q (T) is associated with i? g C J q {T), we obtain: 

LEMMA 5.29: The Janet bundles are associated with R q and we set T r = T(F r ). 

From this lemma we shall deduce the following important but difficult theorem: 

THEOREM 5.30: The first order operators V r : F r _i ->■ F r induce maps V r : T r _i -> T r in 
the deformation sequence: 

— »• Z(6) — »• (7(0) At„At 1 A..At„^0 
which is locally described by finite dimensional vector spaces and linear maps. 

Proof: Any section of F r _i can be lifted to a section A 7 ^ 1 G A r ~ 1 T* (£>J q (T) modulo A r - 1 T*®R q + 
<5(A r_ T* (gi^g+iT* CgT) where the second component is in the image of D. Then we can lift again 
this section to a section A r q ~\ G A r "T*® J q+ i(T), modulo a section of A r "T*®i? 9 +i+L»(A r "T*(g) 
J 9 9 +2 (T)) + A r - 1 T*(g)S' 9+ T*®r where J q q+2 is the kernel of the projection tt|+ 2 : Jq +2 (T) -> J q (T). 
The image by £> r is obtained by applying D : A r ~T* ® J 9+ i(T) -> AT* ® J 9 (T) and projecting 
DA T q +l 6 AT* J q {T) thus obtained to F r while taking into account successively the restriction 
£> : A r_1 T* O R q+1 -> AT* the fact that £> o L> = L> 2 = and the restriction 5 providing 
an element in 5(A r ~ 1 T* ® S q+ \T* ® T). Of course, with such a choice we need to use associations 
with Rq+i even though finally only R q is involved because of the above lemma. The reason is that 
the restriction of D to J" +1 (T) has an image in T* ® J q (T) and woi in T* <g) J q {T). Hence, starting 
with a section of Y r _i C iA_i, we must have: 

Lfe+OA^" 1 G A r -T* ®R q + 5{Ai- 2 T* ® S 9+1 T* T) V^+i G R q+1 = Pl (R q ) 
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Accordingly, looking at the right member in the formula of Proposition 5.13, the section in 
/\ r_1 T* eg) R q is contracted with r — 1 vectors in order to provide a section of R q and the skewsym- 
metrized summation finally produces a section of A r T*®R q . Similarly, as 6(A r ~ 2 T*<g>S q+ iT*(g)T) C 
A r-i T * S(]T * (g, T we should obtain a section of AT* ® ST* T but it is not evident at all that 
such a section is in <5(A*'~T* ® S q+ iT <g> T). For this, let us notice that, if the multi-index / has 
length r — 2 and \fi\ = q, we have the local formula: 

a f i+i i jjdx' 1 A (dx^ A dx 1 ) = —a^+^jjdx^ A (dx 1 A dx 1 ) 
explaining why the following diagram is commutative and exact: 

T* <g> A r - 2 T* ® S q+1 T* { ~^ S A r_1 T* ® S q+ iT* — > 

4. J 4. J 

r* <g> A r - 1 T* <g> SgT* ("IlV 5 A r r®s,r — >0 

where the upper map induced by <5 : A r ~ 2 T* T* — > A r ~T* anc | the lower map induced by 
<5 : A r ~T* T* — > A r T* are both epimorphisms while the left vertical map is induced by 
<5 : A r - 2 T* S q+1 T* -> A r_1 T* ® 5,7*. 
We obtain therefore: 

DL(Z q+2 )A r q -{ - L(£ q+1 )DA r q -l e AT * + ^(A r_1 T* ® S, +1 T* T) 
Now, from the construction of V r already explained, we have: 

F r _! = A r_1 T* J q+ i(T) / (A r ~ 1 T* ® R q+1 + D(A r ~ 2 T* ® J 9 9 +2 (T)) + A^T* ® S, +1 T* ® T) 
It follows that: 

L(^+ 2 )^l € A r - X T* ® R q+1 + D(A r ~ 2 T* ® J g 9 +2 (T)) + A r_1 T* ® S, +1 T* ® T 
and thus: 

DL{i q+2 )A r q -\ e AT* ® i? 9 + 5(A r -T* S g+1 T* ® T ) 

that is finally: 

L(Z q+1 )DA r q -{ e AT* ® i? 9 + 5(A r " 1 T* S g+1 T* ® T ) 
Accordingly, DA r q z\ G AT* ® J g (T) is the representative of a section of T r . 

Q.E.D. 

DEFINITION 5.31: The deformation sequence is not necessarily exact and only depends on 
R q . We can therefore define as usual coboundaries B r (R q ), cocycles Z r (R q ) and cohomology groups 
H r (R q ) = Z r (R q )/B r (R q ) with B r {R q ) C Z r (R q ) C T r for r = 0, 1, ...,n. 

PROPOSITION 5.32: We have B (R q ) = C(Q)/Z(@),Z (R q ) = N(Q)/G and the short exact 
sequence: 

o — ► c(e)/z(e) — > jv(e)/e — > ffo(# g ) — > o 

Proof. The monomorphism on the left is induced by the inclusion C(0) C AT(0) because Z(Q) = 
0nC(0). Then Bo(R q ) is the image of V in To because of Proposition 5.23 and Proposition 5.24. 
Finally, Z (R q ) = N(Q)/Q is just a way to rewrite Corollary 5.27 while taking into account the 
fact that T>\ o V = 0. The cocycle condition just tells that the label transformations induced by 
the normalizcr do not change the structure constants because the Vessiot structure equations are 
invariant under any natural transformation. This result should be compared to Lemma 2.7. 

Q.E.D. 

In order to generalize Proposition 5.26 and to go further on, we need a few more concepts from 
differential geometry. 
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DEFINITION 5.33: A chain £ £' £" of fibcrcd manifolds is said to be a sequence with 
respect to a section f" of E" if im(<&) = kerjn (^), that is with local coordinates (x, y) on £, (x, y') 
on £', (x,y") on £" and y' = $(x, y), y" = *(x, y'), we have 1-(x, $(x, y)) = f"(x), V(x, y) G £. 

Differentiating this identity with respect to y, we obtain: 

— ( X ,^ X ,y))—( X ,y)=0 

and we have therefore a sequence V{£) — > V(£ ) — > V(£ ) of vector bundles pulled back over 
£ by reciprocal images. 

DEFINITION 5.34: A sequence of hbered manifolds is said to be an exact sequence if im($) = 
kerf"(^) and the corresponding vertical sequence of vector bundles is exact. 

PROPOSITION 5.35: If £,£',£" are afEne bundles over A with corresponding model vector 
bundles E,E', E" over A, a sequence of such affine bundles is exact if and only if the corresponding 
sequence of model vector bundles is exact. In that case, there is an exact sequence £ £' — > E" 
which allows to avoid the use of a section of £" while replacing it by the zero section of E". 

Proof. We have successively y' = &(x,y) = A(x)y + B(x),y" = ^(x,y') = C{x)y' + D(x) 
and by composition y" = C{x)A(x)y + C{x)B{x) + D(x) — f"(x). Accordingly, we must have 
C(x)A(x) = 0,C(x)B(x) + D(x) = /"(x),Vx G A and we obtain the following commutative dia- 
gram: 



E ^ E' V J$ E" 

; ; / ; 

£ A £' A £" 
tt | n'l n" |t /" 

A = A A 

If (x, y') e £' is such that C(x)y' + D(x) = f"(x), we obtain by substraction C(x)(y' — B(x)) = 
C(x)(y' — 3>(x,y)) = with (x, y' — $(x, y)) G £?'. Supposing the model sequence exact, we 
may find (x, v) € E such that y' - (A(x)y + B(x)) = A(x)v <^> y' = A(x)(y + v) + B(x) and thus 
(x, y') G im(<£>). The converse is similar and left to the reader. 

Finally, we just need to notice that (x, y") G £" ^ (x, y" - f"(x)) G E" , Vx G A with y" - f"(x) = 
C(x)(y' — B(x)) — C(x)(y' — $(x, y)) where(x,y) G £ is any point over x G A. 

Q.E.D. 

Coming back to the construction of the Vessiot structure equations from the knowledge of the 
generating differential invariants at order q, we have already exhibited the system B\ G J\{T) 
locally defined by affine equations of the form I{u\) = A(u)u x + B(u) = 0. The symbol Hi G 
T* ® V(T) of this system is defined by linear equations of the form A{u)v x = and the last 
proposition provides at once the following commutative and exact diagram of affine bundles and 
model vector bundles. 



0^ Hi — »• T* <g> V{T) — > Ti — >0 

: '■ S II 

0^ Bi — »• Ji(J") — > J"i — >0 

4' \r \" 

T = T = T 

More generally, having in mind the diagram at the end of section 4, we may define the nonlinear 
Janet bundles T r = AT* ® V(J r )/5(A r ~ 1 T* ® Hi) as a family of natural vector bundles over T. 
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The following commutative diagram of reciprocal images generalizes the result of Proposition 5.26: 



section of T r J/f 4-t equivariant section 

X T 

As a byproduct, we may generalize Proposition 5.26 by saying that any infinitesimal change 
of source x = x + t£(x) + ... can be lifted to J> with u T = u T + + ...,v a = 

U a + tt*(x)Mg\1(u)vP + .... 

REMARK 5.36: No classical technique could provide this result because all the known methods 
of computer algebra do construct the Janet sequence " step by step" and never " as a whole" , that 
is from the Spencer operator ([36], p 391). 

PROPOSITION 5.37: The affine bundle: 

A(R° q ) = {x g erg J q {T)/(T* ® R° q + S(S q+1 T* ® T))| Xo - idr G T* ® T} 
is modelled on F\ and associated with R q . 

Proof. First of all, A(R° q ) is modelled on T* ® J°(T)/(T* ® i?° + 5(5 9+ iT* ® T)) = F x and it just 
remains to prove that the affine bundle: 

A q = { Xq G Ci(T)| Xo - i<*r G T* ® T} = id^(Ci{T)) 

is a natural affine bundle over X modeled on the vector bundle C?(T) = T*® J®(T) / 5{S q+1 T* ®T) 

and associated with J q (T), thus with i? g too because T* ® R q is ass ociated with R q . 

For this it is sufficient to observe the transition laws of T* ® J q (T) when x — <p(x), namely: 

t*®t x\Av r = x k A<p l 



T*®J q (T) ^... r „A^-ft./^V r + i- = xjA^ + - + X fe A +lfe ^ 



with replaced by (ii, i q ) when = q. Setting x\ — wc get: 

id T \T* ® Jx(T)) X l ,,A<P r d j <p' = Xj,i9k<p l + dijV 1 

idr\T* ® J,(T)) xk-r.^i^-^^^i^ + = X^<W + ... + 

and thus: 

= [Zq,Xq{v)] + i(v)D£ q +l - Xq(L{Zl)v) 

because Xo(v) — V- If we choose G J q+ \{T) over £ 9 e Jq{T) 1 the difference will be 

i{rj)D(£ q+ \ — Cg+i) = — i(»7)<S(£g+i — Cg+i) — by residue because we are in C\{T). We notice that 
no formal integrability assumption is needed for R q . 

Q.E.D. 

Similarly, we obtain: 

PROPOSITION 5.38: id~ 1 { Ji(II 9 )) is an affine natural bundle of order q + 1 over X, modelled 
on r®J,(T). 

Proof. We only prove the proposition when g = 1 as the remaining of he proof is similar to the 
previous one. Indeed, if x = <p(x), y = ^p{y) are the changes of coordinates on X x Y, we get: 

ni m<p r W = Tf(y)y? 

J x (HO V?,,d i <p r (x)d j <p>(x)+y?d ij <p'-(x) = ^{y)yl J + p§ ¥ {y)y^ ] 

idzHM^i)) y^A^dj^ + dij^ = dk^ylj+dik^y* 
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Then, we just need to set y = x, = Sj and compare to T* ® Ji{T) above. 

Finally, as Ji(II g ) is an affine bundle over H q modelled on T* <g> V"(II 9 ), then id~ 1 ( Ji(H q )) is an 
affine bundle over X, modelled on id- x {T* ® V(U q )) = T* ® id' 1 (V '(II 9 )) = T* ® J q (T). 

Q.E.D. 

PROPOSITION 5.39: We have the following exact sequences of affine bundles over X and 
model vector bundles: 



T* ® R q — > T*®J q (T) 



id-H-MUg)) — ► id- 1 (Mil,)) — > w-^j^t-)) — >o 

■i - "i- 4" 

x = x = x 

Proof. We just need to use 1 <| [i \< q and set y q = id q (x) in order to obtain: 

u r = * T (u,{y),y ll )%u> T {x) 

< = iB^^M + M iyq)yk ^ ^ - L " iu}{x))yh ^ + yk * dkU}T{x) 

in a coherent way with the upper model sequence and the various distinguished sections. 

Q.E.D. 

PROPOSITION 5.40: We have the following exact sequences of affine bundles over X and 
model vector bundles: 

— ► T*®R° — > T*®J°(T) V -^ T*®F — >0 
0^ id T \T* ®R q ) — ► id^{T* ®J q (T)) — ► w-^J^J 7 )) — >0 

4 4' 4" 

X X X 

Proof. As id 9 is a section of then j\(id q ) is a section of Ji(IT 9 ) over id q and thus a section of 
id~ 1 (Ji(II ? )) denoted by If x 9 = (if,Xuj) with 1 <| [i < g is a section of id^(T* ®J q (T)), 

that is a J 9 (T)-connection, it is of course a section of T* <g) J q (T) and we may consider the section 
y q ,i = id q ,i+Xq of id^iJiiUg)). The image in w _1 (Ji(.F)) is u[ = (-L^(cj(x))x^ + fta; T (x)) + 
cW"(x) and we obtain: u[ — diW T (x) <^> ui = ji(w)(a;) x 9 £ id^, (T* ® 

Q.E.D. 

THEOREM 5.41: ~ cj-^Ji^/Bi) 

Proof. The affine sub-bundle £>i C Ji(F) over J 7 is the imager of the first prolongation pi($) : 
— >• Ji(J-) where n g+ i is an affine bundle over II^ modelled on S q +iT* ® V(X x X). 
We may proceed similarly by introducing the affine bundle id^ 1 (n g+ i) which is modelled on 
id^iSg+iT* ® F(X x X)) = 5 9+ iT* ® 1 (V r (-XT x X)) = S q+1 T* ® T and obtain the com- 
mutative and exact diagram of affine bundles and model vector bundles where the model sequence 
is the symbol sequence that has been used in order to introduce F\: 



S q+1 T* 



F 



M-^+i) 
X 



X 



u-\J x {F)/Bi) 

I 
X 



28 



With 1 <| fi \< q, the first morphism of affine bundles is now described by: 

because x does not appear in $ T (y 9 ), in such a way that ul~diU! T (x) = ~^2^^ =q Lj, ,i (ijj(x))y^ +li G 
T* ® F . It just remains to notice that 5 : S q+ iT* ® T — >T* ® S g T* ® T is a monomorphism. 

Q.E.D. 

If x g £ T* ® J 9 (T), with projections xo = G T* ® T and xi G T* ® Ji(T), is a rep- 
resentative of an element c G we may choose a lift Xg+i G T* ® J q+ i(T) and define 
§{x«> Xg} G A 2 T* ® J°(T) by the formula: 

^{Xg,Xj(»7,C) = {Xq+l(v),Xq+l(0} ~ Xq{{Xl{v),Xl(()}) 

PROPOSITION 5.42: The map Xq — > \{Xq,Xq} provides a well defined map c — > j{c,c} 
from elements of A{R° q ) invariant by R° q to F 2 = A 2 T* ® J°(T)/(A 2 T* ® R° q + S(T* ® S q+1 T* ® T)). 

Proof. First of all, if we change the lift Xg+i to Xg+ij then Xg+i — Xg+i G T* Sg+iT* <8> T and 
the difference will therefore be in <5(T* (g) S q+ iT* ® T) with a zero projection in F 2 . 
Then, if we modify Xq by an element M q G 6(S q+1 T* ® T) C T* ® S q T* ® T C T* ® J q {T), 
the difference in J q -\{T) will be only produced by {Xq(v)i Xg(C)} an d wm be {M g (ry), Xg(C)} + 
{Xqiv): M q (()} = 8M q (r], C) = because 5 o 5 = 0. However, when lifting M g at order q + 1, new 
terms will modify the quadratic application and, in particular, even one more if q = 1, namely 
Xi({Mi(0,Xi(»7)} + {Xi(£)> M i (»?)}) = = xi(0) = and such a situation will not 

differ from the general one we now consider. According to the above results, the only pertubating 
term to study is: 

{<rl+l(0,X g +i(v)} + {Xq+iiO^T+M) - M g {{xi(Z), Xi (V)}) G S q T*®T 

and we just need to prove that such a term comes from an element in 8(T* ® S q +\T* ® T). With 
more details and | /i |= q — 1, we get for the 6 = 3 + 3 factors of ^,rp: 

(x r t M +lr+l] - x r t ,M +lr+u ) + (M; +lt+liX ^ - M; +lt+ljX k r .) Ml +lr+lt ( x i 3 - xli) 

Setting: 

Nq+i = {Nl +u+lj . = (x r t M +lr+lj + X r ,M +lr+lt ) x k r ,iMl +u+h ) ef® S q+1 T* ® T 

where the first sum in the parenthesis insures the symmetry in t/j, a tedious but straightforward 
though unexpected calculation proves that all the above terms are just described by 5N q+ i. 
As another proof, we may also consider the 3-cyclic sum obtained by preserving \i but replacing 
(i,j,t) successively by (j,t,i) and (t,i,j) in order to discover that all the terms disappear two by 
two. It is then sufficient to use the exactness of the (^-sequence: 

— > S q+2 T* ® T T* ® S q+ iT* ® T -A A 2 T* ®S q T®T A 3 T* ® S q -iT* ® T 

Hence it just remains to modify Xq by cr° G T*®R q with lift a q+1 G T*®J q+1 (T) and the difference 
will be: 

{Xq+M,cr q+1 (0} + W q+ i(v),X q+ i(C)} + kS(l),^(()] 

-X 9 (KW,Xi(0} + {Xi(r?),<7?(C)})-a 9 ({xi(7?),Xi(C)}) 
As the third and the last terms already belong to R q , we need just consider: 

{Xq+li^^+AOy + Wq+MXq+liOy-XqiSo-^O) 

Now we have: 

i(C)L(a q+1 (v))Xq = L(a° q+ i(ri))Xq(0 Xg^KWK) = K + iM,Xg + i(C)} - X,(i(^W)() 
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and the above remainder becomes: 

i(OL(a q+1 ( V ))x q -i(v)L(a° q+1 (C) £ J°(T) 

because ' 

i(^(i))C - H<ri(Q)v - t°i(v, C) = o 

Finally, L(^)c = 0,V£° £ in Fi 4=> L(^ +1 ) Xg E T* ® R° q + 5(S q+1 T* ® T). Hence, modifying 
<Jq + i if necessary by an element in T* <g) 5 g+ iT* ® T, the above remainder belongs to 
It follows that \{ Xq , Xg} will be modified by an element in A 2 T* ® R° q + 5(T* 5 9+ iT* ® T) that 
will not change the projection in F 2 . 

Q.E.D. 

If we start now with any system R q C J q (T) of infinitesimal Lie equations satisfying [R q , R q ] C 
R q , we may choose any i? g -connection Xq £ id^, (T* <S> R q ) C id^, (T* ® Jq{T)) and obtain by 
projection an element c £ that we use in the following proposition. 

PROPOSITION 5.43: ir^ 1 : R q+1 — ► i? 9 is surjective if and only if L(£ q )c = 0, V£ q E R q . In 
this case, we have \{c, c} = in F 2 . 

Proof. We have Xo(^) = ?7 and we may find £ q+ i E J q+ \{T) such that: 

i(v)L(tq+l)Xq = L(tq+l)Xq(v) ~ Xq(L(tl)v) 

= [£,, X,fa)] + ifaW+i - X^iM £ J 9 °(T) 

As [R q , R q ] C i? 9 and Xq £ r* ® R q , modifying by an element in S q +\T* ® T if necessary, we 
get at once: 

£(&+i)x, e T* ® ^ Lfe)c = £ J 9+1 (T),^ +1 (£ 9+1 ) = C 9 £ i^.I^+i £ T* ® i? 9 

Proceeding backwards in the definition of D given at the beginning of section 4, we obtain therefore: 

L{Qc = £ = p 1 (i 1 , 9 ),7r9+ 1 (e 9+ i) - C 9 £ i? 9 

It follows that 7r| +1 : — > -Rq is surjective and we can choose a lift Xq+i £ 2""* ® over 
X g £ T* ® i? 9 in such a way that ±{x g , Xq] £ A 2 T* ® R° q => |{c, c} = in F 2 . 

Q.E.D. 

REMARK 5.44: As already noticed, the IC L(£ q )c = may be obtained in two steps. First, in 
order to obtain the surjectivity n^ +1 : R q+1 — > R q , we must have L(£ q )c = 0, V£° £ R q , that is c 
wit/si 6e R q -invariant. Second, in order to obtain the surjectivity 7Tq +1 : R q +i — > T, we must have 
L(Xg(0) c — 0, £ T for any i? g -connection Xq- 

The following definition only depends on R q : 

DEFINITION 5.45: A truncated Lie algebra is a couple (R° q ,c) where R° q C J°(T) is such that 
[R° q , R° q ] C i?° and c £ is i?°-invariant with \{c, c] = in F 2 . 

THEOREM 5.46: The infinitesimal deformation of Cart(R q ) = {c £ A{R q v )\c is R q -invariant, 
i{c,c} = 0} is just Z x {R q ). 

Proof. If c t = c + tC + ... is a deformation of c and Xq(t) = Xq + + ■•■is a representative of c t , 
we must have C £ F\ and X° £ T* ® J q (T). Denoting simply by {c,C} = the linearization of 
I{c, c} = 0, we must have { Xq , X°} £ A 2 T* ® #° + <5(T* (g) 5 g+ iT* ® T) and get successively with 
different lifts: 

{x,+i(0,^ifa)}+{*^i(0,x,^ 
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= L( Xq (0)X° q (v) ~ L( Xq ( V ))X°(0 

-i(®DX° +1 { V ) + i(ri)DX° +1 {t) X°({}) Xq () 
= i(,7)L(x,(0)X° - mL(xM)X° q - DX° q+ i(Z,v) X,() 
where we have used the relation: 

i(v)Hx g (0)X° g = L( Xq (Z))X° q (v) X° q (L( Xl (Ov) 
with i(xi(C)) 7 ? = [£; v] + i(v)DXi(0 an d the relation: 

i(77)£>x 9 ° +1 (0 - iCO^V^) + x ? ° ([£,??]) = (£,»?) 

that is finally: 

{x g , r?) + ^°+i(^) = i(l)V^,° - i(OV„^ - x ,() 

Now L(£ q )C = 0,V^ g Rg = R° q © x 9 (T) =► V C C = in F U V£ € T V € X° e T* ® i?° + 
<5(S , 9+ iT* © T) and thus: { Xqi X^} + £>X° +1 g A 2 T* ® Modifying if necessary the lift 
we need only DX° +1 g A 2 T* ® a result leading to L(l q )C = 0,V^ g that is C g Ti and 
V 2 C = in F 2 . 

Q.E.D. 

As the F r only depend on i?° both with the sub-bundles E r C F r , we finally obtain ([37], p 721): 

COROLLARY 5.47: The deformation cohomology H r {R q ) only depends on the truncated Lie 
algebra (R q ,c). 

It finally remains to study a last but delicate problem, namely to compare the deformation 
cohomology of R q to that of R q+ \. First of all, as R q is supposed to be involutive in order to 
construct the corresponding Janet sequence with Janet bundles F r , then R q +i is of course involu- 
tive too and we may construct the corresponding Janet sequence with Janet bundles F' r . The link 
between the two Janet sequences is described by the next theorem. 

THEOREM 5.48: There is the following commutative diagram with exact columns and locally 
exact top row: 
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Proof. We may use the following commutative and exact diagram in order to construct successively 
the epimorphisms \& r : Fj. — >■ F r with kernels if r by using an induction on r starting with r = 0. 























i 




1 




1 


1 


-> 


.9<?+r+l 


-> 


S g+r+ iT* ® T 


-> 




-> ir r ->■ o 




1 




1 




1 


1 


-> 


R q +r+l 




-/ 9 +r+i(r) 


-> 




-> f; -> o 




I 




1 7 r q+r+1 




1 Jr(*o) 




-> 


R q +r 




Jq+r(T) 


-> 


J r (F ) -»-... 


-> .f; -> o 




I 




I 




1 


1 





















31 



It follows that the Janet sequence for T>' projects onto the Janet sequence for T> and we just 
need to prove that the kernel sequence, made up by first order operators only, is exact. For this, 
introducing the short exact sequences : 

-> g q+ r+i -+ S q+r+1 T* <g> T ->• h r+1 -> 

with h r+ i C 5 r+ iT*(g)Fo the r-prolongation of the symbol hi C T*®F of the system B\ C Ji(F ), 
image of the first prolongation J q+ i(T) — > J\{Fq) of the epimorphism J q (T) — > F used in order to 
define P. Here hi is identified with the reciprocal image of the symbol of B\ C JiiJ 7 ) by to. Using 
the Spencer operator D and its various extensions, we obtain the following commutative diagram 
where it is known that the vertical D-sequences are locally exact ([36], [49]). 
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Now the left vertical sequence is exact because, applying the <5-sequence to the previous short 
exact sequence, then hi is involutive whenever g q is involutive. The exactness of the top row finally 
follows from a diagonal chase in the last diagram because the other rows are exact by induction 
according to the previous diagram. 

Q.E.D. 

In order to use the last theorem in a natural way, we shall, for simplicity, restrict to the use of 
the classical Lie derivative by using the formulas: 

= iMO,v q ] v£ g t, y Vq e j q {T) 

aOivM = WOvM + {v g , £(0C q } 

DC(0 = £(0D 

which are direct consequences of the fact that the algebraic bracket, the differential bracket and D 
only contain natural operations. Accordingly, when restricting to an involutive system R q C J q (T) 
of infinitesimal Lie equations, we have to preserve the section oj of the natural bundle T in order 
to construct the Janet sequence and we obtain at once: 

PROPOSITION 5.49: We have the formulas: 

C(£)V = VC(Q, C(Z)V r =V r C(Q, V£ee,Vr = 0,...,n. 

EXAMPLE 5.50: Whenever V is a Lie operator already defined by the Lie derivative of oj with 
respect to a vector field, we get: 

(C(0v - vc(0)v = c(QC(v)«> - m v\ = - ^) V = C(v)£{t)u = £(v)V£ = o 

Such a result may be applied at once to all the known structures. The case of a symplectic structure 
with n = 2p, oj € A 2 T* , det(oj) ^ 0, dco = is particularly simple because a part of the correspond- 
ing Janet sequence is made by a part of the Poincare sequence and it is well known that the Lie 
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derivative £(£) commutes with the exterior derivative d (See [37], p 682 for more details). 

The two following propositions will be obtained from the diagram of the last theorem by means 
of unusual chases in the tridimensional diagram obtained by applying £(£) for £ G 9 to the diagram 
of the last theorem. However, before providing them, it is essential to notice that a short exact 
sequence — > F — > F' — > F" — >• of bundles associated with R q may only provide in general the 
exact sequence — > T — > T' — > T" by letting £(£) acting on each bundle whenever £ G 9, unless 
one can split this sequence by using a natural map from F" to F' . For example, in the case of the 
short exact sequence — > J q (T) — > J 9 (T) — > T — > one can use j q : T — > J q (T). We also recall 
that an isomorphism Z'/B' ~ H' ^ H = Z/ B induced by maps Z' — >• Z and B' ^ B does not nec- 
essarily imply any property of these maps which may be neither monomorphisms nor epimorphisms. 

PROPOSITION 5.51: There is an isomorphism = H (R q+1 ) ~ H (R q ) = H Q . 

Proof. We shall provide two different proofs: 

1) First of all, this result is a direct consequence of the short exact sequence provided by Proposi- 
tion 5.32. 

2) Let us prove that there is a monomorphism — > H' — > i?o- First, if a' 6 Tq C Fq is killed 
by V[ and such that * (a') = a = Vrj for a certain »j G T(T) C T, then a = * (a') = *o(^'»?) 
and thus V (a' - V'rj) = 0, that is a' - V'r] = b' G K with V[b' = V[a' - V'{D'r) = 0, that is 
b' = => a' = V'r]. 

Let us then prove that there is an epimorphism H' —> H n — > 0. For this, if a G T G F is killed 
by Vi, we may find a' G Fq such that a = ^o(a') and £(£) a ' 6 ^o- It follows that c' = 2?^a' G i^i 
with £>2 C ' — 0) that is c' = 25^6' for a certain &' G i^o- Accordingly, modifying a' if necessary 
by replacing it by a' — b' , we may suppose that V[a' = with ^o(a') = a and we just need 
to prove that C(£)a' = 0. Indeed, we have * (£(O a ') = £(£) a = => £(£)a' £ K o and thus 
V[£(^)a' = £(^)V[a' = => £(£) a ' = because the restriction of X>i to K is a monomorphism. 

Q.E.D. 

PROPOSOTION 5.52: If R q C J 9 (T) is involutivc, the normalizer N(Q) of 9 in T is defined 
by the involutive system R q +i G J q+ i(T) with involutive symbol = g q+ i G 5 g+ iT* ® T, 

defined by the purely algebraic condition {R q+ \, R q +i} G R q . Moreover, Z (R q +i) — Z (R q ) — 

Rq+l/Rq+l- 

Proof. With Z = Z n (R q ) and Zq = Z (R q+ i) as before, let us prove that there is an isomor- 
phism ->■ Zq ->■ Z Q ->■ 0. For this, if a' G Z„ is such that * (a') = 0, then a' G K with 
V[a' = a' = because the restriction of T>[ to K is injective and we get a monomorphism 
-> -> Z -> 0. 

Similarly, if a G Zo, using the previous proposition, we can find a' G Fq with 2?^ a' = and 
*o(a') = a. It follows that £(£)a' = d' e K with X>i<f = £>i£(£K = C{£)V[a' = d' = 
£(£) a ' = =>■ a' G and we get an epimorphism Zq — >• Z — >• 0, that is ~ Z n (R q ). 

Now we have C(£ q+ i)r] q = {£ q+ i,r) q+ i} + i(£)Dri q+ i G i? 9 for any lift T] q+ i G J q +i{T) of 77, G J 9 (T). 
However, cocycles in Z are also killed by V\ and the representative r\ q G must be such that 

Dr/ q+ i G T* (g) + (5(S' 9+ iT* ® T). Therefore , modifying r^+i if necessary without changing rj q , 
we may suppose that Dr/ q+ i G T* a result leading to the condition {£ q +i, G i? g => 

r/ 9 +i G -Rg+i- As we already know that = fif g +i, introducing the projection i? g C J q {T) of 
.Rq+i G J q+1 (T), we obtain Z Q (R q ) ~ ^/i?, ~ R q+1 /R q+1 and similarly Z (i? 9+ i) ~ R q+2 /R q +2 
in the following commutative and exact diagram: 
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Chasing in this diagram or counting the dimensions in order to get: 

dim(R q+ 2) - dim(R q+ i) = dim{R q+ 2) - dim(R q+ i) = dim(g q+2 ) = dim{g q+ 2) 

it follows that 71"^ ■ Rq+2 — > Rq+iis an cpimorphism. 
Next, applying the formula of Proposition 5.6, we get: 

{£ q +i,i{0Dv q +2} = i(()D{Z q +2,V q +2} - {i{0D£ q+ 2,ri q +i} ^R q ^ Dr lq+2 G T* <g> R q +i 

It follows that R q +2 = Pi(R q +i) because both systems project onto R q +i and have the same sym- 
bol g q+2 = g q+2 = Pi(g q +i) = Pi(g q +i)- 

Finally, as g q +i — g q +i is an involutive symbol when g q is involutive and Tr q +1 '■ R q +2 — > Rq+i is 
surjective, then R q +i is involutive because of the Janet/Goldschmidt/Spencer criterion of formal 
integrability ( [23] , [36] , [39] , [49] ) . 

Q.E.D. 

PROPOSITION 5.53: There is a monomorphism ->• H[ ->• H 1 . Accordingly, Hi(R q+1 ) = 
whenever Hi(R q ) = and thus R q +i is rigid whenever R q is rigid and we may say that is rigid. 

Proof. If a' G C F[ is killed by V 2 and such that ^i(a') = a = T>\c for a certain c G T C F , we 
may find d G Fq such that * (c') = c and thus a = ^i(a') = 2?i* c' = ^iD[d =4> *i(a'-X>ic') = 
^> a' - P^c' = 6' G As before, we get X^6' = V' 2 a' - V 2 V' lC ' = => 3d' G K with 

D^<f = 6' and thus a' = T>T>\{d + d'). Hence, modifying c' if necessary, we may find d G Fq such 
that ^o(c') = c with a' = V[d and it only remains to prove that d G Tq. However, we have 
C(£)d G K a because £(£)c = and V[C(^)d = £(£)£>! c' = £(£)a' = 0. As the restriction of V[ 
to Kq is injective, we finally obtain £(£,)d = as we wanted. 

Q.E.D. 

Despite many attempts we have not been able to find additional results, in particular to com- 
pare H2{R q ) and H2{R q +\). However, integrating the Vessiot structure equations with structure 
constants Ct instead of c, we may exhibit a deformation u>t of uj as a new section of T . As T> 
depends on ji(cj), the V r will also only depend on uj and various jets. Replacing uj by uj u we shall 
obtain operators T>(t) and T> r (t) for < r < n. Therefore we may apply all the techniques and 
results of section 2 with only slight changes. 



6 EXPLICIT COMPUTATIONS: 

EXAMPLE 6.1: (f .1 revisited) This is a good example for understanding that the Janet bundles 
are only defined up to an isomorphism, contrary to the natural bundles. In this example, we have 
of course the basic natural bundle T = T* ®V = T*x X -.-X X T* (n-times) with V = W l . This is 
a vector bundle (even a tensor bundle) and we may therefore identity F with Tq — V(F) and Fq. 
However, in this case, gi = R\ = and we get F Q = J$(T)/R% = T* <g> T. Similarly, we should 
g etF r = A r T*®T*®T/(5(A r " 1 T*®52T*®T) = /\ r+1 T* ®T because we have the exact sequence : 
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A r_1 T* O S^T* -A A r T* (g>T* ^ A r+1 T* — > 
Using geometric objects, we should obtain: 

and thus fi = A 2 T* <g> V. On the contrary, using the solved form £ — > (V£)* = <9i£ fe + 
£f a k (x)d r ul(x) providing at once a unique i?i-connexion, a covariant derivative V and thus a 
well defined V-sequence leading to F r = A r+1 T* ® T and thus Fi = A 2 T* ® T. In fact, i/ie expk- 
nation is a confusion between the Janet sequence and the Spencer sequence because we have indeed 
C = R 1 =T^ CV+i ^ F r ~ A r+1 T* ® T. 

The sections of AT* ® V invariant by Ri are of the form A T T± (Jr uj c ' 1 A ... A u ar while the sections 
of A r T* ® T invariant by R\ are of the form a k A T (Ji (Jr ijJ c ' 1 A ... Auj (7r and are thus different though 
depending on the same constants A T a a , providing therefore a unique cohomology. From the 
Maurer-Cartan equations duj T = —c T pa u p Aw" where we change the sign of the structure constants 
fo convenience, we obtain terms of the form Ac which are exactly the first terms in the definition 
of the Chevalley-Eilenberg operator. Now, looking at the terms containing a, we obtain (care to 
the minus sign): 

(d ia k T + a r T a k d r wf)A T dx i = a r T a k (d r uj? - d i uj a r )A T dx l = -a k c T pa A p uj a 
and recover the second terms of this operator. 

It is however essential to notice that the two approaches are totally different. In particular, ele- 
ments in A(R®) are 1-forms with value in some vector bundles and it is therefore a pure chance 
that they could become 2- forms in this example (See electromagnetism and gravitation in [43]). 
We may also compare the concept of "change of basis" of the Lie algebra with the "label transfor- 
mations". Indeed, using the Lie form, we get (care again to the minus sign): 

Hi ~ to a k d r uij — a k d r ujf = —oj?d r a k 

a p d r u)J + u) T k d r a k = 
«• r (a*a>r) = O 
<=> ujJ = a T a u°~ , a = est 

Accordingly, the sections of To C T* ® V are of the form A t g uj° with A = est and it is rather 
extraordinary that two such different approaches can provide the same result. 
Finally, we study 6, Z(<d),C(0) and N(Q) in this framework. First of all, is described by a 
linear combination with constant coefficients of the infinitesimal generators {0 T = O l T di}. Then 
we must have — £%if + £, r d r n k = 0, whenever £f + ^ r a k T d r (jj\ = 0, that is a k d r ulrf + d r r\ k = 
or d r (ujjr] 1 ) = and thus rf — X 7 a\ with A = est {reciprocal distribution). We already know 
that [apjOv] = c pa a a with our choice of sign. The isomorphism between G and C(6) is there- 
fore a pure chance. As for N(Q), it follows from Corollary 5.27 that the defining equations are 
(£(£)w)[ ee u)l{x)di^ r + £, r d r uil(x) = A T a ujf t {x) where A — est is a derivation of the Lie algebra Q 
with the same structure constants c. 

EXAMPLE 6.2: (1.2 revisited) Let us study the equivalence relation Hi ~ co by asking first that 
R\ = R\. We must have u> r j£l + &ir£,j = <S=> to r j£l + w ir£J = 0. Setting £ij = ui r jS,i, it amounts 
to check that ui r jUj rt £,i t t + uii r uj rt (;j.t = £ij + = and we must therefore have: 

uj rj uj rt 8l + Gj lr ijj rt b' s j - cD rj ui rs 8 t i - uj lr uj rs Sj = 

Contracting in s and j, we get nui r jOJ — ui rs oj rs Sj = 0, that is ui r jUj rt = a{x)8 t j w t3 = a(x)oJij. 
Substituting in = we finally get ui = auj,a — est and the group of label transformations is 
just the multiplicative group. It then follows that 7 = 7 but such a property can be checked directly 
from the solved form of the second order equations £(£)7 = 0. It follows that p — p and thus 
c = (l/»)c and, in any case, N(&) is defined by the infinitesimal Lie equations £(0 W = Auj with 
cA = as there is no Jacobi condition on the single structure constant c. Accordingly, N(Q) = 6 
if c ^ and dim(N(@)/Q) = 1 if c = 0. This is the reason for which the Weyl group is the normal- 
izer of the Poincare group, obtained by adding the generator x l 9, of the dilatation on space-time. 
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It is important to notice that the Galilee group is of codimension 2 in its normalizer obtained 
by dilatating separately space and time, another reason for which the Poincare group cannot be 
obtained from the Galilee group by a continuous deformation. As for C(O), using again we 
have to solve rfu^&j = <=>■ rf&j = 0, V&j + = and thus C(6) = => Z(Q) = 0. 

EXAMPLE 6.3: (1.3 revisited) This example is by far the most difficult to treat in dimension 
n = 2p + 1 ([37], p 684). In the present case when n = 3, this is one of the best examples where the 
Lie equations obtained by eliminating p, namely £3 — x 3 £f = 0, £2 — x3 £,i + x3 £i — ( x3 ) 2 £i — £ 3 = 
could be used without even knowing about the underlying geometric object. However, in this 
case R\ is not involutive and we must start afresh with the involutive system C R\ for con- 
structing the canonical Janet sequence. We first notice that there is only one CC for the new 
involutive system i?i C J\(T) of infinitesimal Lie equations and thus surely no Jacobi condition 
for the only structure constant c. It thus remains to study the inclusions Z(Q) C 6 C N(Q). 
Using Rl C Ji(T) for C(9), we have to solve r/ s ^ — whenever w r £[ - ^lo^ — and 
thus i] s cu s ^ = => r/ s uj s =0. As w ^ 0, changing coordinates if necessary, we may sup- 
pose that uii 0. It follows from the involutive assumption that at least one jet coordinate 
of each class is parametric and in fact, changing coordinates if necessary, we have in all the ex- 
amples presented pri = {?|, ^3,^2} =^ P ar — {£1 > £1 > £1 > £l> £f> £3}- Hence, choosing successively 
as unique non-zero parametric jet, we obtain n 1 = 0, ?? 2 = 0, 77 s = => r\ = and thus 
C(6) = Z(©) = 0. As for N(@), we have to study the equivalence u> ~ w, that is to study 
when we have w r £[ — = <^> w r £[ — 5 w i£r- Though it looks like to be a simple algebraic 

problem, one needs an explicit computation or computer algebra and we prefer to use another more 
powerful technique ([37], p 688). Introducing the completely skewsymmetrical symbol e = (e' 1 ' 2 * 3 ) 
where e tlt2 ' 13 — 1 if (i-ii^is) is an even permutation of (123) or —1 if it is an odd permutation and 
otherwise, let us introduce the skewsymmetrical 2-contravariant density ui 1 ^ — e n u>k- Then one 
can rewrite the lie equations as : 

Rl - {x)Cr ~ - \^ + Cdr^ 3 (X) = 

and we may exhibit a section = ut ls A rs with A rs — A sr and thus = 0. It is important to notice 
that det(uj) = when n = 2p+l, contrary to the Riemann or symplectic case and w cannot therefore 
be used in order to raise or lower indices. As we must have Ri — the same section must satisfy 
(u, r iuj is +uj ir u> js )A rs = 0, VA rs = A sr , and we must have (w r J ' uj is +tj"'uj js ) + (u)^ ' uo lr +Q ls uj^ r ) = 0. 
Setting s = j, we get w^'a/- 7 ' = Co^uj'^ (D 4 - 7 (x) = a(x)u) v (x). Substituting and substracting, we 
get a/ 3 ' (x)^ r d r a(x) = => a(x) = a = est or similarly uji(x) = a(x)u>i(x) =>• u>i(x)£, r d r a(x) = 
a(x) = a = est because ui ^ and one of the components at least must be nonzero. It 
follows at once that one has N(T) = {/ 6 aut(X) | ji(/) _1 (w) = aui,a 2 c = c}. Accordingly, 
N(Q) — {£ G T C(^)lo = Auj, Ac = 0} and © is of codimension 1 in its normalizer if c = or 
N(Q) = © if c ^ 0. For example, in the case of a contact structure with c = 1, we have N(Q) = 
but, when ui = (1,0,0) => c = 0, we have to eliminate the constant A among the equations 
c^C 3 + d2xi 2 — diS, 1 = —2A, d^ 1 = 0, 82s, 1 = and we may add the infinitesimal generator x l di of 
a dilatation providing A = — |. 

EXAMPLE 6.4: (1.4 revisited) This is by far the most interesting example. First of all, contrary 
to all the previous examples, we have at once T = J" = T*x x A 2 T*, T\ = A 2 T*x x A 3 T*,J 7 2 = 
A 3 T* and we may identify F r with T r for r = 0, 1, 2, 3 in order to obtain the Janet sequence: 

with dim(T) = 3,dim(F ) = 6,rfim(iq) = 4,dim(F 2 ) = 1 and 3-6 + 4-1 = 0. Using local 
coordinates (u 1 ,?/ 2 ,^ 3 ) for T* and (u 4 ,u 5 ,u 6 ) for A 2 T* in F , we obtain W — [W\ = u 1 -^- + 

u2 7&+ u3 l&> W 2= ui mF+ u5 (&+ u6 (&}- Withw = W). wehavePC- (£(0", >C(0/3)"and 
recall the iG made by the Vessiot structure equations da = c'/3, d/3 = c"a A (3 with 7=aA/3^0. 
As before, we may easily find the label transformations {a — aa, f3 = b/3) =>• (c! — %c', c" = -c") 
and obtain T = (Aa, B/3) C F a with A, B = est. Similarly, we get Ti = (C/3, C"j) C F 1 with 
C, C" = est and the induced map V 1 : T -> Ti : (A, B) -> (c'(A - B) = C", -c"A = C") in a 
coherent way with the linearization of the Vessiot structure equations: 
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Ada - c'B/3 = c'(A - B)/3, Bd/3 - c"Aa A /J — c"Ba A /? = -c"A 7 

We obtain therefore A(6) = {£ G T | £(£)a = Aa, £(0/8 = B0, c'(A - B) = 0, c",4 = 0} and thus: 
c = (0, 0) => dim(N(Q)/e) = 2, c = (1, 0) A = B rfim(Af(6)/6) = 1, c = (0, 1) A = 
dim(N(Q)/Q) = 1. The only Jacobi condition c'c" = provides c"C + c'C" = in any case. For 
example, in the unimodular contact case, we must add x 1 d\ + x 2 <9 2 to 9 in order to obtain N(Q). 
As for the centralizer C(0), we must look first for r\ G T such that rf{£ — 0, G i??. For 
this, multiplying /3 r j£l + Pir£j = by rji and contracting on i, we get {r) l Pi r )£j — when- 
ever a r ^l = 0. Accordingly, we obtain rf fii r — L(x)a r as 1-forms. Now, as n — 3, we may 
introduce the pseudo-vector (J3 1 = /?23,/? 2 = /?3i,/3 3 = /812) transforming like a vector up to 
a division by the Jacobian A and thus, using the volume form 7 = a A /3, it follows that 
/3 = (/3 1 /7, /3 2 /7, /? 3 / 7 ) is a true vector field. As /3i r /3 r = because n = 3, we obtain by contraction 
L{x)a r (5 r = L(x)(a A ,8/7) = L(ar) = => n l fi vr = -> r/ fe = K{x)(3 k . But a,/3,7 and thus /3 are 
invariant by any (66 and thus i\T(x) = K = est, that is C(9) = {0 G T|?] = if /3, if = est}. In 
all the three special sections considered, we have simply C(9) = {{K, 0, 0) G T|K = est}. 

We obtain finally the folowing recapitulating picture: 



(K) — ► — ► (C'.C") — ► (-D) — ►() 

(if/3) — ► (Aa,B/3) — ► (C"/3,C" 7 ) — ► (D7) — ► 

— > 9 — ► T T* x x A 2 T* A 2 T* x x A 3 T* A 3 T* — > 

£ — ► ($,*) — > (J7,V) — > (W) — >0 

The purely differential lower part is describing the operators involved in the Janet sequence, 
namely and successively: 

v £— > (£(£)<* = £(00 = *) 

2?! ($, *) — ► (d$ - c'* = U,cM- c"(/3 A $ + a A *) = V) 

2> 2 (?/, v) — >• (d?7 + cT = W) 

The purely algebraic upper part is induced by these operators acting on the second line which 
is describing the invariant sections of the respective Janet bundles and we obtain the linear maps: 

V (K) — ■» (0, c"K) 

X>i (A, B) — ► (c'(A — B) = C, —c"A = C") 

V 2 (C, C") — ► (c"C + c'C" = D) 



only depending on the structure constants c = (c', c"). 
7 CONCLUSION: 

The work of E. Vessiot, motivated by the application of the local theory of Lie pseudogroups ([51]) 
to the differential Galois theory, namely the Galois theory for systems of partial differential equa- 
tions ([52]), has been deliberately ignored by E. Cartan and followers ([8], [26]). As a byproduct, 
the Vessiot structure equations, introduced as early as in 1903, are still unknown today. Similarly 
and twenty years later but for other reasons related to his work on general relativity ([6]), in par- 
ticular his correspondence with A. Einstein, Cartan did not acknowledge the work of M. Janet on 
the formal theory of systems of partial differential equations (Compare [7] and [23]). Accordingly, 
the combination of natural bundles and geometric objects with differential sequences, in particular 
the linear and nonlinear Janet sequences has been presented for the first time in a rather self- 
contained manner through this paper. The main idea is to induce from the linear Janet sequence, 
considered as a linear differential sequence with vector bundles and linear differential operators, a 
purely algebraic deformation sequence with finite dimensional vector spaces and linear maps which 
may not be exact. We have shown that the equivalence problem for structures on manifolds has 
to do with the local exactness of the Janet sequence at F while the deformation problem of the 
corresponding algebraic structures has to do with the exactness of the deformation sequence at 
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Ti = Y(i<i), the space of invariant sections of F\, that is one step further on in the sequence. This 
result explains why the many tentatives done in order to link the deformation of algebraic struc- 
tures like Lie algebras with the deformation of geometric structures on manifolds have not been 
successful. Meanwhile, we have emphasized the part that could be played by computer algebra in 
any effective computation and hope to have opened a new field of research for the future. Finally, 
ending with a wink, the reader must not forget that one of the first aplications of computer age- 
bra in 1970 has been done in the deformation theory of Lie algebras where a few counterexamples 
could only be found in dimension greater than 10, that is with more than 500 structure constants !. 
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